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^ ■ Abstract: 



We analyze a completely integrable two-dimensional quantum-mechanical model that emerged 
in the recent studies of the compound gluonic states in multi-color QCD at high energy. The 
5^ \ model represents a generalization of the well-known homogenous Heisenberg spin magnet to 
infinite-dimensional representations of the SL{2, C) group and can be reformulated within the 
quantum inverse scattering method. Solving the Yang-Baxter equation, we obtain the i?— matrix 
for the SL{2, C) representations of the principal series and discuss its properties. We explicitly 
construct the Baxter Q— operator for this model and show how it can be used to determine the 
energy spectrum. We apply Sklyanin's method of the Separated Variables to obtain an integral 
representation for the eigenf unctions of the Hamiltonian. We demonstrate that the language of 
Feynman diagrams supplemented with the method of uniqueness provide a powerful technique 
for analyzing the properties of the model. 
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1. Introduction 

In this paper, we study a completely integrable quantum-meclianical model of A^ interacting 
spinning particles in two-dimensional space. The model can be thought of as a generalization 
of the well-known periodic one-dimensional XXX Heinsenberg spin chain magnet ll|, §, ^ - 
each particle defines the site of the spin chain and the interaction occurs between spins of two 
neighboring particles. The spin operators associated with each particle are the generators of an 
infinite-dimensional principal series representation of the SL{2, C) group and that is why we shall 
refer to the model as a noncompact spin magnet. 

The motivation for studying such models comes from two different and seemingly unrelated 
physical problems. The first of them has to do with the Regge phenomena in QCD. As was 
shown in 0, ^ , the high-energy behavior of the scattering amplitudes in the generalized leading- 
logarithmic approximation is governed by the contribution of the color-singlet compound states, 
built from an arbitrary number {N = 2,3,...) of gluons and defined as ground states of the 
effective QCD Hamiltonian. The effective Hamiltonian acts on the two-dimensional transverse 
coordinates of the gluons and exhibits remarkable properties of integrability @, |^. Namely, 
in the limit of the multi-color QCD, the wave function of the A^— gluon compound state turns 
out to be identical to the ground state of noncompact Heisenberg spin chain model with the 
number of sites equal to the number of gluons, A^. The quantum space in each site of this model 
is parameterized by two-dimensional gluon transverse coordinates and corresponds to the spin 
s = principal series representation of the SL{2, C) group. At A^ = 2 the ground state is known 
as the BFKL Pomeron ||^. At A^ = 3 it gives rise to the Odderon and for higher A^ the 
corresponding ground states define the unitarity corrections to the scattering amplitudes. In 
spite of a lot of efforts, the solution of the Schrodinger equation for A^ > 3 gluon compound 
states still represents an outstanding theoretical problem in QCD. 

Another motivation for studying the noncompact spin magnets comes from the problem of 
solving the quantum completely integrable models defined on an infinite-dimensional phase space. 
The best known example of such models is given by the one-dimensional Toda chain model |^. 
One of their main features is that they can not be solved using the conventional Algebraic 
Bethe Ansatz (ABA) method [|l], H, ^ and one has to rely on more advanced methods |1T0| , [11 



The noncompact two-dimensional SL{2, C) spin magnets belong to the same class of models 
and their solution represents a theoretical challenge. As we will show in this paper, these two- 
dimensional models have many features in common with the known one-dimensional integrable 
models. Namely, they admit the same i?— matrix representation as XXX Heisenberg compact 
spin chain and part of their spectrum can be reconstructed using the ABA approach. At the 
same time, similar to the Toda model, the quantum space of the system is infinite-dimensional 
and, in general, the eigenstates do not admit the highest weight representation. This means that 
the ABA approach can not provide the full set of the eigenstates of the model and, therefore, it 
is not complete. 

Our approach to solving the Schrodinger equation for the SL{2, C) spin magnet relies on the 



methods of the Baxter Q— operator |]T0[ and the Separation of Variables (SoV) developed by 
Sklyanin [|n|]. The former allows to determine the energy spectrum of the model, while the latter 
provides an integral representation of the corresponding eigenfunctions. In recent years, both 



methods have been developed and successfully applied to solving the Toda chain model [TT, |T2|, 



and high spin generalizations of the homogenous XXX Heisenberg SL{2, M.) spin magnet [Q. In 
this paper, we shall extend these results to the SL{2, C) Heisenberg spin magnets. 



The central point of our analysis is the representation of different operators in the model 
(the i?— matrix, Hamiltonian, Baxter Q— operator etc.) by their integral kernels defined on the 
two-dimensional plane. In this way, remarkable integrability properties of the model that are usu- 
ally expressed in the form of operator identities, like Yang-Baxter equations, are translated into 
(complicated) integral relations between the corresponding kernels. Adopting the language of the 
two-dimensional Feynman integrals, one can associate the kernels with particular Feynman dia- 
grams and develop the diagrammatical representation for the above identities. One of the main 
findings of this paper is that, using the diagrammatical approach, one can establish integrability 
properties of the SL(2, C) spin magnet without doing any actual calculations by applying two 
elementary identities between the Feynman diagrams known in QCD as the "uniqueness rela- 
tions" [T^, r^|. Following this approach, we find the solution of the Yang-Baxter equation for the 



i?— matrix corresponding to the SL(2, C) principal series representation, construct the Baxter 
Q— operator and the unitary transformation to the Separated Variables and, finally, establish 
different relations between them that allow to solve the model. 

The paper is organized as follows. In Section 2 we introduce the Hamiltonian of the model and 
show that it admits the i?— matrix representation. We obtain the expression for the i?— matrix 
by solving the Yang-Baxter equation for the principal series SL{2, C) representation and discuss 
its general properties. This allows to prove a complete integrability of the model and establish 
its symmetry properties. Section 3 is devoted to the construction of the Baxter Q— operator for 
the SL{2, C) spin magnet. Our approach is similar to the one used before for the Toda chain 



[^ and the homogenous SL{2,M.) Heisenberg spin magnets \T^]. Examining the properties of 
the obtained expressions, we establish the relation between the transfer matrices of the model 
and the product of the Baxter Q— operators. This relation allows to express the Hamiltonian of 
the model in terms the Q— operator and reduce the original Schrodinger equation to the problem 
of finding the solutions to the Baxter equation on the eigenvalues of the Q— operator under the 
additional conditions imposed by the analytical properties of the Q— operator and its asymptotic 
behavior at infinity. In Section 5, we apply Sklyanin's method of Separation of Variables to 
obtain an integral representation for the eigenfunctions of the Hamiltonian. We construct the 
unitary transformation to the separated variables and demonstrate its relation with the Baxter 
Q— operators. We obtain the quantization conditions on the separated variables and calculate 
the integration measure on the space of the eigenstates in the SoV representation. Concluding 
remarks are presented in Section 6. Appendix A provides a detailed description of the method 
of uniqueness which is intensively used throughout the paper. In Appendix B, we describe the 
properties of the SL{2, C) transfer matrices including the fusion identities and their relation 
with the Baxter Q— operators. The relation between the obtained integral representation of the 
eigenstates and the Algebraic Bethe Ansatz method is discussed in the Appendix C 

2. The quantum noncompact spin chain model 

2.1. Definition of the model 

The noncompact spin chain model is a quantum mechanical system of A^ interacting particles 
on a two-dimensional (x,?/)— plane. In high-energy QCD, this plane corresponds to transverse 
gluonic degrees of freedom. It becomes convenient to define the position of the particles on the 



plane by introducing complex holomoprhic and antiholomorphic coordinates {k = 1, ■■■,N) 

Zk = Xk + iyk, Zk = Xk-iyk- (2.1) 



We associate with each particle a pair of mutually commuting holomorphic and antiholomor- 
phic spin operators, Sa and Sa , satisfying the standard commutation relations [Sa ,S^] = 

ieafi^5^"'S\ and similarly for Sa ■ They act on the quantum space of the fc— th particle, V^^**'"**), 
and can be represented as the following differential operators 

'S'o = Zkdz^ + Sfc , S- = -(9^^ , 5*+ ^ = zldzt, + 2skZk , 

S'i^ = Zkd-z, + Sfe , S. = -d,, , Sf = zld,, + 27skZk , (2.2) 

so that (5('=))2 = {Sf^f + {Sfs^'^^ + Si'^^Sf V2 = Sk{sk - 1) and similar for the antiholomorphic 
Casimir operator [S^'^^'Y- The spin operators defined in this way are the generators of the unitary 
principal series representation of the SL{2, C) group 

^{Zk, Zk) -^ ^'{zk, Zk) = {czk + d)~^'^{czk + dy'^^^izl 4) , (2.3) 

where 

_ azk + h _ azk + b 

czk + a czk + a 

with k = 1, ...,N, ad — he = ad — he = 1 and the complex parameters Sk and Sk specified below 
(see Eq. (|2.13| )). In what follows we shall assume that the model is homogenous and the particles 
have the same spin, Sk = s and Sk = s ioi k = 1, ..., N . 

The Hamiltonian of the model, TiN-, describes the interaction between A^ noncompact SL{2, C) 
spins attached to the particles and it has the following general form 

Hn = Hn + Hj^ , [Hn, if at] = , (2.5) 

where H^ and H^q act on the holomorphic and antiholomorphic coordinates of the particles, 
respectively, and therefore commute. The (anti)holomorphic Hamiltonian is given by the sum of 
two-particle Hamiltonians describing the nearest neighbor interaction between the corresponding 
(anti)holomorphic spins with periodic boundary conditions 

N N 

Hn = Y, H{Jk,k+i) , H^ = Y, H{Jk,k+i) , H{J) = i,{l -J)+ i,{J) - 2V^(1) (2.6) 
fc=i fc=i 

with %Ij{x) = dliaT{x)/dx being the Euler function and Jn,n+i = Jn,i- Here, Jk,k+i and Jk,k+i 
are defined through the Casimir operators for the sum of the spins [| 

Jk,k+i{Jk,k+i - 1) = (5^'^ + 5^'+'))' (2.7) 

with Sa = Sa , and Jk,k+i is defined similarly. 

This paper is devoted to solving the Schrodinger equation 

Hn ^(2*1, Z2, ..., zn) = En ^(ii, ^2, •••, zn) (2.8) 



3 



Although thi s equation defines the operators Jk,k+i up to substitution Jk,k+i ^ 1 — Jk.k+i, the two-particle 



Hamiltonian (2.6) is invariant under this transformation. 



with the eigenstates \l/(zi, ..., zat) being single-valued functions on the plane z = {z,z), normal- 
izable with respect to the SL{2, C) invariant scalar product 

||^||2= I d^Zid^Z2...d^ZN\'^{zi,Z2,...,ZN)\^ (2.9) 

with d'^z = dxdy = dzdz/2. The equation ( p.8| ) has appeared for the first time in the analysis of 
the high-energy asymptotics of the partial waves in the multi-color QCD 0. It was found [^, 0], 
that it possesses the set of mutually commuting conserved charges whose number is large enough 
for the Schrodinger equation ( p.8|) to be completely integrable. 

The total spin of the iV— particle system is one of the conserved charges. Indeed, the Hamil- 
tonian (|2.6|) is a function of two-particle Casimir operators and, therefore, it commutes with the 
operators of the total spin Sa = ^^ Sa and Sa = ^^ Sa , acting on the quantum space of the 
system V^ = V^(*i''^i) y^'^^''^^^ ® ... V^^^''^'^\ This implies that the eigenstates can be clas- 
sified according to the irreducible representations of the SL(2,C) group, V^'^''^\ parameterized 
by the spins (h, h). The eigenstates "ifiz), belonging to V^'^'^^ are labelled by the center-of-mass 
coordinate zq and they can be chosen to have the following SL{2,C) transformation properties 

N 

^(4 - z',) = {czo + df^czo + df"- \[{czu + df'-{czk + df'^^{zu - zo) (2.10) 

fc=l 



with Zq and zq transformed in the same way as Zk and Zk, Eq. ( |2.4| ). As a consequence, they 
diagonalize the Casimir operators corresponding to the total spin of the system, S'^ = Sq + 
{S+S^ + S_5+)/2, 

(52 _ /^(/^ _ 1)) ^^^^^ ^^^ _ _^ ^^) ^ (^2 _ iQ^ _ ^)^) ^(^^^ ^^^ __^ f^) = . (2.11) 



The complex parameters (s^, Sk) and (/i, h) entering ( |2.2| ) and ( |2.1CI| ) parameterize the irre- 



ducible S'L(2,C) representations. For the principal series representation they satisfy the condi- 
tions g 

Sk-Sk = n,, , Sk + {skY = 1 (2.12) 

and have the following form 

Sk = — ^ + tJys, , Sk = — ^ + tus, (2.13) 

with Usi^ being real and rig^, being integer or half-integer.Q The spins {h, h) are given by sim- 
ilar expressions with n^^ and Usf. replaced by Uh and Uh, respectively. The parameter Ugf. 
has the meaning of the two-dimensional Lorentz spin of the particle, whereas z/^^, defines its 
scaling dimension. To see this, one performs a 27?— rotation of the particle on the plane, 
z -^ zexp{2'n'i) and z -^ zexp{—27Ti) and finds from ( [^.10 ) that the wave function acquires a 



phase \l/(zfc, Zk) — > (— l)^"^'fc\I'(2;fc, Zk). For half-integer Ug^ it changes the sign and the correspond- 
ing unitary representation is spinor. Similarly, to define the scaling dimension, s + s = 1 + 2ii>Sf., 

^ Throughout the paper we indicate by bar symbol the variables belonging to the antiholoniorphic sector. 
Notice that, in general, the corresponding variables in two sectors are not complex conjugated to each other. 
^ Since the unitary representations labeled by the spins (s, s) and ( 1 — s,\ — s) are unitary equival ent a nd are 



related to each other through the intertwining relation (see Eq. ( 3.76| ) below) one can choose Us^. in ( 2.13 ) to be 
nonnegative. 



one performs the transformation z ^ \z and z — > \z. We recall that for the homogenous spin 
chain one takes Sk = s and Sk = s for arbitrary k. In what follows, for the sake of simplicity, we 
will not consider the spinor SL{2, C) representations and choose n^^ in ( |2.13| ) to be nonnegative 
integer. 



We notice that the holomorphic and antiholomorphic spin generators ( |2.2| ) as well as the 
Casimir operators (|2.7|) are conjugated to each other with respect to the scalar product (p.9| ) (see 
footnote to dg^)) 

[Sa \ = —Sa , [Jfc,fc+l] = 1 — Jk,k+l ■ (2-14) 

This implies that H]^ = Hj^f in Eq. ( |2.5| ) and, as a consequence, the Hamiltonian (|2.5| ) is hermitian 
on the space of functions endowed with the scalar product (|2.9|), TY]^ = Hn- Thus, the energy 



E]s[ in ( p.8|) is real and the corresponding eigenstates are orthogonal to each other with respect to 
( p.9|) . Since the Euler ■?/'— function has poles at negative integer arguments, the holomorphic and 
antiholomorphic Hamiltonians, Eq. ( |2.6| ), are unbounded operators. One can verify however [^ 



that, thanks to the properties of the principal S'L(2,C) series, Eqs. ( p. 121) and ( p. 131) , the poles 



are cancelled in their sum ( p.5|) leading to the Hamiltonian TYtv which is bounded from below. 

2.2. i?— matrix 

Let us show that the model defined in the previous section can be described using the i?— matrix 
approach |jl| . For this we will need the expression for the i?— matrix acting on the tensor product 
of two SL{2,£,) representations, V"*-*^'*^^ ® y(s2,s2) ^ Th.e general expression for the i?— matrix 
is well known for the Heisenberg spin magnets in the case of arbitrary higher spin 5'L(2,M) 



representations [^, |j and to the best of our knowledge its SL{2, C) analog has not been discussed 
in the literature. 

To find the i?— matrix for the infinite-dimensional SL{2, C) representations, we introduce the 
Lax operators in the holomorphic and antiholomorphic sectors 

u + iSo iS^ 



L.(u) = u + i(a-S)-, _^^ ^ _ .<,^ 

im = S + .(..S)=(^+£^« ,;L%) (2,15) 

with u and u being arbitary complex parameters and a^ being Pauli matrices. These operators 
act on the quantum space ^(^'*) and coincide with similar expressions for the Lax operator in 
the Heisenberg spin chain model |jl], |], |^. We define the i?— matrix by requiring that the Lax 
operators have to satisfy the commutation relations 

Ls^{v)Ls.^{v + ^)-R(si,si),(s2,s2)(«, u) = R(si,si),(s2,s2){u, u)Ls^{v + u)Ls-,{v) . (2.16) 

Here, the operator -R(si,si),(s2,s2)('"; '^) depends on two spectral parameters, u and u, and it acts 
on the tensor product ^(■^i-^i) V^(^2,S2) labeled by the SL{2, C) spins (si, si) and (s2, S2). 



2.2.1. Yang-Baxter equation 



Solving tlie Yang-Baxter equation ( |2.16|) , we use tlie metliod proposed in |]19|, [2^]. It is based on 



the representation of the i?— operator by its integral kernel |j 

[RisuSi),{s2,S2){u,u)'^] (^1,^2) = / (i^Wi / (fw2Ru,u{^l,Z2\Wi,td2)'^{Wi,W2), (2.17) 

where Zi = {zi,Zi) and Wi = {wi,Wi) are two-dimensional vectors with the (anti)holoniorphic 
coordinates defined in (|2.1|) and \l/(wi,W2) is an arbitrary test function belonging to ^(*i'*i) (g) 
Yis2,s2)_ ^Jq f[T^(^ ^]^Q kernel, Rufi{zi, Z2 \wi, W2), we require that the i?— matrix defined in this way 
has to satisfy the Yang-Baxter equations ( |2.16| ) and, in addition, Ru,u{zi,Z2 \wi,W2) has to be 
a single- valued function on the two-dimensional plane. As we shall see in a moment, the latter 
condition imposes constraints on the possible values of the spectral parameters u and u. 

The dependence of the kernel of the /2— matrix on the holomorphic coordinates is fixed by 
the first relation in ( ^2.161 ). Applying its both sides to the same test function "^{^1,^2), one 



substitutes the Lax operators and the i?— operator by their expressions, Eqs. (|2.15|) and (|2.17|) , 
and integrates by parts in the r.h.s. using the identity 

[R{u,u)S^^^'i/] (ZI,Z2) = / rfV / d^W2RuA^l,^2 \wi,W2) {Si''\wi)'^(wi,W2)) 

= - [d^wi [d^W2{S^^-'^\wi)R^4zi,Z2\wi,W2))^iw,,W2i2.18) 



Here, Sa^ (wi) and Sa * (t/^i) denote the differential operators acting on the holomorphic coor- 



dinates wi. They and given by Eqs. (|2.2|) with the spin Sk replaced by si and 1 — si, respectively. 
In this way, the first relation in ( |2.16| ) can be replaced by the following matrix equation 

{v + ia- S^''\zi)) {v + u + ia- S^'^\z2)) RuA^u ^2 \wi, W2) 
= {u + v-ta- S^^'''\w2)) {v-ia- S^^-''\w^))Ru,u{zi, Z2 \w^, W2) , (2.19) 

which leads to an overcompleted system of the differential equations on the kernel Ru,u- It turns 



out that the system has a unique solution. The simplest way to find it is pO| to project the both 



sides of (|2.19| ) by the vector (— ^2, 1) from the left and by the vector (1, Z2) from the right and. 



then, match the coefficients in front of powers of u and v. One obtains two first-order differential 
equations on Rufi, which fix the dependence of the kernel on the holomorphic coordinates Zi and 
Wi up to arbitrary prefactor depending on the specral parameters. Repeating similar analysis for 
the second equation in ( p.l6| ), we restore the dependence of the kernel on the antiholomorphic 
coordinates. The resulting expression for the solution to ( p.l6| ) is given by the product of the 
holomorphic and antiholomorphic kernels 

Ru,u{^l,Z2\Wl,W2) = Pji{u,u) 

X{W2 - Z^y^-''+''-\zi - Z2)-'''-''-''+\wi - W2y'^''^''-\Z2 - 101^"+''-''-' (2.20) 
X(W2 - Zi)*"-'"^+'"^"Vzi - Z2)-'''-''-''^\wi - W2)-'''^''^''-Uz2 - u;i)™+«i-«2-l 



^Here, we assume that there exists the region of the spectral parameters, u and u, in which the integral is 
convergent. For u and u outside this region the kernel of the i?— operator is defined by the analytical continuation. 




-1 + S^ + $2 



+ lU 



(Z-W) 



Z-W) 



Figure 1: Feyninan diagram representation of the kernel of the i?— matrix, Eq. (|2.2(]| ). Each side 
of the square represents a two-dimensional propagator and the arrow indicates in which order 
the difference of coordinates is taken. 



with the normalization factor p^ {u, u) being an arbitrary function of the spectral parameters. 
The obtained expression for the kernel of the i?— matrix can be translated into the language 
of two-dimensional Feynman diagrams as shown in Fig. |I[ As we will demonstrate below, this 
diagrammatical representation becomes extremely useful in solving the model. 

Examining the expression for the kernel, Eq. ( p^.2UD , we notice that Ru^^Zi, Z2 |wi, W2) acquires 
a nonzero monodromy as Zi encircles the points Z2 and W2 on the plane and similarly for other 
arguments of the /?— matrix. For the kernel to be a single-valued function, the corresponding 
monodromies should cancel in the r.h.s. of ( |2.20| ). In general, for the function of the form 
{w — zi)"'{w — zi)"- this condition amounts to the Lorentz spin a — a to be integer. Applying 
the same condition to ( |2.20|) and taking into account the expressions for the spins si and S2, 
Eqs. (|2.13| ) with n^^ and ris^ integer, we find that the spectral parameters u and u have to satisfy 
the additional condition 

(2.21) 



i[u 



u 



n 



with n being an integer. 

Using the expression for the kernel (|2.20 ) it becomes straightforward to show that the 
i?— matrix defined by Eqs. ( p.l7[ ) and ( |2.20| ) satisfies the Yang-Baxter equation 



i?12(u, u)Ri3{v, V)R23{V -U,V ~u) = R23{v -U,V - u)Ri3{v, v)Ri2{u, u) 



(2.22) 



where we used a shorthand notation for Rkn{u, u) = -R(sfc,Sfc),(s„,s„) (w, u). Going over to the integral 
representation ( |2.17|) , the product of the i?— operators in the l.h.s. of ( p.22|) has the following 
kernel 



Ruiu, u)Ri3{v, v)R23{v -u,v -u) 



[Zi,Z2,Z3\Wi,W2,W3) 



(f'yi <f'y2 d^y3 RuA^i, Z2\yi, y2)Rv,v{yi, Z3\wi, y3)Rv-u,v-u{y2, y3\w2, W3) . (2.23) 



One finds similar expression for the r.h.s. of ( p.22|) . The proof of ( |2.22|) can be carried out 
diagrammatically as shown in Fig. ^ without doing any calculations. Replacing each i?— matrix 
by a square (see Fig. 0) one represents the l.h.s. of ( 2.22|) by the left Feynman diagram in the 
lower line in Fig. |^. This diagram consists of three squares and each pair of squares shares 
a common vertex. The central triangle formed by the sides of three squares turns out to be 
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3 X 





Figure 2: Diagrammatical proof of the Yang-Baxter relation based on the "uniqueness" relations. 
Integration over the position of internal vertices (fat points) is tacitly assumed. 



"unique" (see Appendix A for the definitions) and it can be transformed into the "star" diagram, 
using the uniqueness relation (|A.8|) . Then, three new stars with the centers located in the vertices 
of the central triangle are also unique and can be transformed back into triangles giving rise to 
a hexagon diagram. Performing the same set of transformations on the r.h.s. of ( p.22| ) (the left 
diagram in the upper line in Fig. ^, one arrives at the same hexagon diagram, thus proving the 
Yang-Baxter equation. 

Following similar procedure one can show that the i?— operator (|2.20|) satisfies the T— inversion 
relation (see Fig. |) 



^(si,si),(s2,S2)\'^y '^j-^(si,si),(s2,S2)( 



vr p^[u,u)pj,[-u,-u) 



a[si — S2 + iu) a{S2 — si + iu) aisi — S2 — iu) a[s2 — si 



(2.24) 
X 1, 



lU] 



where the function a(x) is defined in ([A.3|) and 1 stands for the identity operator on V^'^^''^^' 
y{s2,s2) _ Defining the normalization factor p^{u^u) we fix, for later convenience, the coefficient 
in front of the identity operator in Eq. (|2.24|) to be equal to unity 



1 



Pif{u, u) = — a{si — S2 + iu) a{s2 



Sl 



lU] 



(2.25) 



2.2.2. Eigenvalues of the i?— matrix 



It follows from the Yang-Baxter equation ( p.l6| ) that the i?— matrix commutes with the sum the 
SL{2, C) spins and, therefore, it is invariant under the SL{2, C) transformations (|2.3| ) and (|2.4| )n 

[Ru{u, u), S'^' + S'^'] = [Ruiu, u), SL'' + SL''] = . (2.26) 

To see this, one takes the hmit in ( 2.16| ) as w — > oo for u = fixed and similarly for the antiholomorphic sector. 



yi 



^1 yi 



= 5(2^- y^) S(Z2- y2) 



y2 



y2 ^2 y2 



Figure 3: Diagrammatical proof of Eq. (|2.24| ). Two lines connecting wi and W2 cancel each other 
and the integration over the position of these points is performed using the chain relation ( [A.7| ). 






Figure 4: Diagrammatical calculation of the eigenvalues of the i?— matrix. 

This implies that the i?— operator is a function of the Casimir operators J12 and J12 defined 
in ( p.7|) . In order to find the explicit form of this function or, equivalently, to calculate the 
eigenvalues of the i?— matrix, one has to decompose the tensor product V^^^^i'^i) V"(^2,S2) ^y^j- ^j^g 
irreducible components V^'^''^^ and take into account that the operators, J12, J12 and Ri2{u,u) 
are diagonal on V^'^'^' simultaneously. 

The projector onto the principal series spin— (/i,/i) representation, n'^'^''*\ is defined as 



\^ih,h)l 



£w, d^W'^U^'^'^Uw^ - 



Wi — Z,W2 — Z 



|^(S1,S1),(S2,S2)| 



Wi,W2 



(2.27) 



where \l>('^''^)(i') belongs to V^^''^^ and the kernel ]l^'^''^\wi — z,W2 — z) is given by |T7 



n^'^'^^Vu;! - 



wi -z,W2-z) = (wi - ^2)'^+'^+'^"' {W2 - zY^-^-'' {z - wiY'-''-'' 

X (Wi - W2f+''+'^-^ (W2 - z)*"2-ft-Si /^ _ ^ xs-i-h-s-2 



(2.28) 



Requiring this kernel to be a well-defined function on the plane, one finds that the possible values 
of the spins {h, h) are given by general SL{2, C) expressions ( p. 13 ) 



n = h tUh ■ 



n = h tUh 



(2.29) 



with Uh arbitrary real and Uh integer. One can verify that the projectors Il^^''^\zi — zq, Z2 — Zq) 
are orthogonal to each other with respect to the scalar product ( ^.91 ) for different Zq and the 
SL{2, C) spins [h, h) such that n^ > 0. 
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By the definition, the projector diagonalizes the Casimir operator (5*1 + 6*2)^ = h{h — 1) and 
{Si + 82)^ = h{h - 1) leading to 

n^^'^) J12 = h n^'*''^) , n^'*''^) J^ = h u^^'~^^ . (2.30) 

To find the eigenvalue of the i?— matrix we calculate the product of the operators U^'^''^^ Ri2{u, u). 
Replacing the operators by their kernels, Eqs. ( |2.28| ) and ( p.20|) , one gets 



d''zid^Z2 n('^'^)(zi -Z,Z2- z)Ru,uiZi, Z2\Wi, W2) = Rh hiu, u) U^''~^\wi - Z, W2 - z) , (2.31) 



where the r.h.s. is fixed up to factor Rf^j^{u,u) by the SL{2,C) transformation properties 
The calculation of the integral in the l.h.s. can be perform diagrammatically as shown in Fig. |^. 
The corresponding diagram is obtained by gluing together the triangle (n*^^'''^) and the square 
{Ru,u) along the line connecting the points zi and Z2. In the resulting diagram, these points are 
the centers of two stars, which turn out to be unique. Subsequently applying the uniqueness "star- 
triangle" relations, one obtains the triangle diagram, which is equal to the projector multiplied 
by a c— valued factor Ry^j^{u,u) depending on the 5'L(2,C) spins h and h. Replacing them by 
the corresponding operators, Eq. (|2.3CI| ), one obtains the operator form of the i?— matrix 

_ V{s2 — si + iu)V{l + si — S2 + iu) Til — h — iu)V(h — iu) 
Rhrh{u,u) = — . .^,, , r^ X 



r(s2 — si — m)r(l + Si — S2 — iu) r(l — h + iu) T{h + iu) 
Ri2{u,u) = RhjXu,u) _ _ . (2.32) 

h=Ji2,h=Ji2 

Remarkably enough, the obtained expression for the /?— matrix is factorized into the product 
of the holomorphic and antiholomorphic operators. Each of them formally coincides with the 
well-known expression for the i?— matrix for high spin representations of the SL{2,M.) group ||18|| . 
However, the important difference with the latter case is that the spectral parameters in two 
sectors, u and u, are not arbitrary anymore and have to satisfy the additional condition ( p.21| ). 



2.2.3. Unitary i?— matrix 

For the homogenous spin chain, Si = S2 = s and Si = S2 = s, the general expression for the 
i?— matrix, Eq. ( |2.32 ), simplifies to 



T{iu)T{l + iu) T{l-Ji2-iu)T{Ji2-iu) 
Ri2[u,u) = — — . , , r^ X — — ■ x-p/ T — r— T- [Z.66) 

r(— ■j'u)r(i — lu) r(i — J12 + m) r(Ji2 + m) 

This operator acts on the tensor product V 0V with V = \/(*'*) and has the following properties. 
At u = u = the /?— matrix (|2.33|) coincides with the permutation operator P12I 

Ru{0, 0) = (-1)1+-^--?- = -p,2 , (2.34) 

which is defined on K ® V^ as 

Pl2^(^l,i2) = ^(i2,^l). (2.35) 



^ Since the spectral parameters have to satisfy the condition ( |2.2l| ), the hmit of ( ^.33 ) at m = u = has to 
calculated by putting u = u and sending u — » afterwards. 
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To verify ( p.34| ), one projects its both sides onto II^'^'^) and takes into account ( p.29|) 



Taking into account (|2.14|) , one finds that the i?— operator satisfies the following relations 

[Ri2{u,u)]^ = Ru{-u\-u*), Ri2{u,u)Ri2{-u,-u) = l. (2.36) 



According to ( 2.36| ), the i?— matrix is not a unitary operator on V ® V for arbitrary u and u 
satisfying (|2.21|) . However, there exists a region of the spectral parameters 



u = u 



(2.37) 



in which the unitarity holds 

[Ri2{u,u)f Ri2{u,u) = l. (2.38) 

Combining ( p.3?l ) together with ( |2.21|) , we find that the i?— matrix is a unitary operator on V"® V" 
for the spectral parameters of the general form 



tn 



m 



u 



u 



(2.39) 



with V real and n integer. 



2.3. Complete integrability 



Let us show that the i?— matrix (|2.33| ) defines a completely integrable quantum-mechanical sys- 
tem with the Hamiltonian given by Eqs. (|2.5| )- (|2.6| ). Applying the i?— matrix approach and fol- 



lowing the standard procedure [0, ^, H, we define the family of the transfer matrices T 
parameterized by the spins (sq, Sq) and acting on the quantum space of the system V 



(«o,so), 

N 



U, U) 



T 



(so, so), 

N ' 



U,U) 



-'-^(scso) [^{so,So),(si,Si){'Uy '^)^iso,So),{s2,S2)\'^J ''^J • • •-^(so,So),(s]v,Sjv) l"^) ^ 



(2.40) 



Here, the trace is taken over the auxiliary SL(2,C) representation space V^(*0'*o). We recall that 
the spin chain is homogenous, so that si = ... = sn = s and si = ... = sn = s. The spins (sq, sq) 
have the form ( p. 131 ) and, in general, they are different from the spins {s,s). Substituting the 
i?— matrices in (|2.40| ) by their integral representation, ( |2.17 ) and (|2.2CI|) , one can evaluate the 



-i{S0,S0) 



kernel of the transfer matrix T)J°'*° (m,m) as A^— fold convolution of the i?— kernel with periodic 
boundary conditions. The diagrammatical representation of the transfer matrix T 



(so, so), 

N 



U,U) IS 



shown in Fig. |^. 

Invoking the standard arguments [^ |^, H, one finds from the Yang-Baxter equation, Eq. (|2.22|) , 
that the transfer matrices form the family of mutually commuting SL{2, C) invariant operators 



[T^^'^^''^\u,u),T^;;^'''^\v,v) 



[Sa,T^^"''^i 



u,u) 



[So.,T^^'"'^'\u,u)] 







(2.41) 



and, therefore, serve as the generating functions of the integrals of motion and the Hamiltonian 
of the model. The latter is obtained from the fundamental transfer matrix T^'^'^\u, u), for which 
the auxiliary space V^^'^^ coincides with the quantum space of a single particle, V, and the 
i?— matrices entering ( 2.401) are given by ( 2.33 ). 

Examining the expansion of the fundamental transfer matrix T^'^''^\u,u) around the origin, 
u = u = 0, and using the properties of the i?— operators, Eqs. ( |2.33|) and ( p.34| ), one finds 



n 



N 



d 

du 



- ln-j(s,s-), 



^N 



U, U) 



H.2 + ... + H, 



u=0 



N-1,N 



+ H, 



NA 



(2.42) 
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Figure 5: Diagrammatical representation of the transfer matrix T^°''^° (u,u) defined in (|2.4U|). 



Tlie left- and rightmost vertices are located at the same point and integration over its position 
as well as over the position of the remaining internal vertices is implied. 



where the two-particle Hamiltonian H^ k+i is given by 



H 



12 



d 

UU u=0 

V^( J12) + ^{l - J12) - 27/^(1) + V^( J12) + ^{l - J12) - 2V^(1) 



(2.43) 

with the Casimir operators J12 and Ju defined in ( |2.7] ). Here, we substituted the i?— matrix by 
its operator expression (|2.33|) and used (|2.34|) . Comparing ( p.42| ) with ( p^SD we conclude that 
two Hamiltonians are identical. 

To identify the total set of the integrals of motion of the model, one constructs the auxiliary 
holomorphic monodromy matrix |I], |^, |^ 



Tn(u) = Li(u)L2(u)...Ln(u) 



An{u) Bn{u) 
Cn{u) Dn{u) 



(2.44) 



and similarly for antiholomorphic monodromy operator T]sr{u). Replacing the Lax operators by 
their expressions L^ = Ls^{u), Eq. ( |2.15| ), one obtains the operators A^, ..., -Dat in the form of 
polynomials in the spectral parameter. Their asymptotics at large u is given by 









(2.45) 



^Af 



where Sa = ^k=i^°' ^^ ^^^ total SL{2,C) spin of N particles. These operators act on the 
quantum space of the system and they are related to their antiholomorphic counterparts as 

[Am{u)]^ = Am{u*) , [Bn{u)]^ = Bm{u*) (2.46) 

and similarly for the remaining operators Cn and D^. 

Taking the trace of the monodromy matrix ( |2.44D we define the auxiliary transfer matrix 



,A^ 



,N-2 



tN{u) = An{u) + Dn{u) = 2u'' + q2u''-^ + ... + g^v 
and similarly for ii\f{u). Combining together (|2.46| ) and (|2.47|) one finds that 



[tN{u)]^ = iN{u*) 



ii = ik 



(2.47) 



(2.48) 



13 



for k = 2, ..., A^. Here, the operators qk {%) are given by certain linear combinations of the product 
of k spin operators. They can be rewritten using ( p. 2]) as A;— th order differential operators acting 
on (anti)holoniorphic coordinates. For instance, 

N 

q2 = -Yl 2(^^'^ ^^"^) = -S' + ^^(^ - 1) = E(^^ - ^nf^^'^^d^.d^^izk - Zr,f' + 2N s{s - 1) , 

k>n k>n 

(2.49) 
where S"^ = h{h — 1) is the Casimir operator corresponding to the total spin of A^ particles, 
Eq. (^J^. 



It follows from the Yang-Baxter equations ( p.l6|) and (|2.22|) that the auxiliary transfer ma- 



trices t]\f{u) and tj\f{u) commute with the transfer matrices T^^° ''^°\v , v) and, in addition, satisfy 
the same relations ( p.41|) . Together with Eqs. ( 2.47 ) and ( p. 42 ) these relations imply that the 



operators {qk, qk) for k = 2, ...,N form the set of 2N — 2 mutually commuting SL{2, C) invariant 
integrals of motion 

[Hn, qk] = [Hn, qk] = [Hn, S^] = [Hn, S^] = (2.50) 

and 

[qk, qn] = [qk, qn] = [qk, S^] = [qk, Sa]=0. (2.51) 

Two additional operators can be added to this set due to the SL{2, C) invariance of the Hamil- 



tonian ( p.50|) . It is convenient to choose them as particular projections of the total spins Sa and 

N N 

P = iS^ = ~iY,d.,, p = iS^ = -iY,d,,, (2.52) 

fc=i fc=i 

which have the meaning of the total momenta of A^ particles. Then, the eigenstates of the 
Hamiltonian with a definite value of the momenta p = {p, p) are given by 

^p.,{,,5}(£i, Z2, ...,zn)= j d^zo e^^»P+*^"°P ^(fi - zo, z^ - zo, ...,zn - zo) (2.53) 

with io being the center-of-mass of the system. 

Thus, the Schrodinger equation ( p. 81) possesses the set of 2A^ mutually commuting conserved 
charges, pand {qk,qk} {2 < k < N), and, therefore, is completely integrable. This implies that, 
firstly, the Hamiltonian of the model can be expressed as a function of the integrals of motion 

Hn = 'HN{q2,q2;--;qN,qN) (2.54) 

and, secondly, the wave function \E'p-{g g}(zi, Z2, ■■■, z^) can be defined as a simultaneous eigenstate 
of the integrals of motion. The corresponding energy levels -E{g,q} can be obtained from ( p.54| ) by 
replacing the operators {q^, qk] by their corresponding eigenvalues. Notice that the Hamiltonian 
can not depend on the momentum operator p due to the SL{2, C) invariance ( |2.50| ). The explicit 
form of the dependence (|2.54|) will be established in the Section 3.5 using the method of the 
Baxter Q— operator. 

2.3.1. Special case: A^ = 2 

At A^ = 2 the Schrodinger equation ( |2.8| ) can be solved exactly. In this case, the Hamiltonian of 
the model is given by the two-particle kernel ( p.6|) , H2 = 2H{Ji2) + 2H{Ji2) and its spectrum 
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can be found by diagonalizing simultaneously the Casimir operators Ju and Ju, as well as the 
momentum operators p and p. Taking into account Eqs. (|2.28|) and ( p.30[) and putting Si = S2 = s 
one finds the eigenstates as 



^p,ihM^l,Z2) = j ^2^0 6^^°^+^^"°^ ^,:d^l - Zo,Z2 - Zo) , 

^^,j,(fi - Zo, Z2 - Zo) = [zi - Z2f~'"'[zi - zo\-\z2 " Zq]"'' , (2.55) 

where the notation was introduced for \zj — za,-]" = {zj — Z]^'^{zj — Zk)°'. Here, the spins h = 
(1 + nh)/2 + iuh and h = {1 — nh)/2 + iuh are the eigenvalues of the SL{2, C) Casimir operators 
Ji2 and Ji2, Eq. ( p.ll|) . One verifies that ( p.55|) satisfies ( p.lO|) . The corresponding energy is 
equal to 



EN=2{h, /i) = 4 Re [iIj{1 - h) + tp{h) - 2^p{l)] = 8 Re 



^(^^J^ + zz.,)-^(l) 



(2.56) 



Minimizing this expression with respect to integer n^ and real u^ one finds that the ground state 
corresponds to /i = /i = 1/2, or equivalently rih = Vh = ^ 

Ef^^ = -lQ\n2. (2.57) 

For s = and s = 1 the expressions ( |2.55D and ( |2.5ti| ) are well known in QCD as defining 
the spectrum of the two-gluon color singlet compound states. The ground state energy (|2.57| ) 
determines the intercept of the BFKL Pomeron [Q. 

2.4. Discrete symmetries 

The Hamiltonian ( |2.42| ) is invariant under the cyclic and mirror permutations of the particles. 
The generators of these transformations, P and M, respectively, are defined as follows 

[P*](fl,...,fAr_l,i'Ar) = ^(^2, ..., ^7V, ^l) , [M ^] (fi, ..., 2^-1, ^Af) = ^(^W, ^Af-1, ..., ^l) 

(2.58) 
so that [H7v,P] = [UnM] = 0. Obviously, the operators P and M are identical at iV = 2. As 
we will show below, this symmetry allows to establish some general properties of the spectrum 
of the model. 

By the definition (|2.58| ), the operators P and M do not commute and satisfy the following 



relations 

P^ = M2 = ]1, pt = p-i = p^-i, M"f = M, PM = MP"^ = MP^~^ (2.59) 

or equivalently (MP)^ = 1. As a consequence, the operators P and M can not be diagonal- 
ized simultaneously and, therefore, the eigenvalues of the Hamiltonian could possess, in general, 
definite quantum numbers only with respect to one of them. 

Let us examine the action of the permutations on the integrals of motion, qk and qk, or 
equivalently on the auxiliary transfer matrices tN{u) and iisfiu). Using the definition ( |2.44| ) and 
( p.47|) , we find that the transfer matrices are invariant under the cyclic permutations, FHn{u)F = 
t^iu), while under the mirror permutation they are transformed as 

Mtjv(«)M = tT{LN{u)...L2{u)Li{u)) 

= {-If ti {Li{-u)L2{-u)...Ln{-u)) = {-iftNi-u), (2.60) 
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where the second relation follows from the property of the transposed Lax operator, Eq. (|2.15|) , 
L*(u) = —a2L{—u)a2. Replacing the auxiliary transfer matrix by its expression ( |2.47D , we find 



Mgfc = (-1)VM, Pgfc = gfcP (2.61) 

and similar relations hold for the antiholomorphic charges. Since the Hamiltonian is invariant un- 
der the mirror permutations, it has to satisfy the following relation as a function of the conserved 
charges, Eq. ( p.54[ ), 

Hiqk, Qk) = MHiqk, qk)M = 7^(MgfcM, Mg^M) = n{i-lfqk, (-l)'gfc) • (2.62) 

This implies that, firstly, the eigenstates of the Hamiltonian corresponding to two different sets 
of the quantum numbers, {g^, Qk} and {{—l)''qk, {—^)''<jk}, have the same energy 

ENiQk, Qk) = En ((-l)'gfe, {-Ifqk) (2.63) 

and, secondly, all energy levels of the Hamiltonian except those with g2fe+i = ^2fc+i = (/c = 
1,2, ...) are (at least) double degenerate and the corresponding wave functions are related as 

^P,{(-l)^gfc,(-l)*<?fe}(^l' ^2, •••, Zn) = [M%,{g,,g-j] (Zi, Z2, ..., Zn) = %,te,g-,} (^7V, •-, ?2, ?l) • (2.64) 

One concludes from ( ^.61| ), that among two operators, P and M, only the first one commutes 
simultaneously with the Hamiltonian and the conserved charges, and therefore, it is diagonalized 
by the eigenstates \E'p'{gg}. The corresponding eigenvalues define the quasimomentum 6 of the 

state 

[P%{,,,-}](zi,...,^jv) =e^'^^'^") %,,,-}(zi,...,^^). (2.65) 

The complete integrability of the model implies that the quasimomentum is a function of the 
total set of the integrals of motion {q, q}. Since P^ = 1, its eigenvalues are quantized as 

e{q, q) = 27r^ , for fc = 0, 1, ..., iV - 1 . (2.66) 



Applying the operator of mirror permutations, M, to the both sides of ( |2.65| ) and taking into 
account (|2.64| ) and ( p.59| ), we find that the quasimomentum of the eigenstate \E'p>{(_i)fcqj. (_i)fcqj.} 
is equal to —6{q, q). This leads to the following relation 

e{{-lfqu,{-lfqk) = - e{qk,qk). (2.67) 

As a consequence, the eigenstate of the Hamiltonian with the quantum numbers q2k+i = q2k+i = 
(k = 1,2,...) has a vanishing quasimomentum and, therefore, is symmetric under the cyclic 
permutations of A^ particles. 

Using the solutions to (|2.65|) , one can construct the eigenstates of the operator of mirror 
permutations M 

^^$^ -, = ^^^ ^pi, o| , M ¥^^ _, = ± ¥^^ _, . (2.68) 

p{q,q} 2 Pi^'^T ' p{g,<?} p{q,q} ^ ' 

Although these states do not diagonalize the integrals of motion, {g, g}, they are the eigenstates 
of the Hamiltonian with the same energy EN{q,q)- We find from (|2.59| ) and (|2.65| ) that the 
operator of cyclic permutations P acts on them as 

P 4tL-} = I (^^' ± ^~'' M) %^.?} ■ (2.69) 
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Notice that at 6* = and 9 = n the states \I/L j --. diagonahze the operator P and, therefore, have 
a definite parity with respect to the cychc and mirror permutations simultaneously. 

Defining the eigenstates of the model, one has to choose between two (equivalent) sets of the 
states, ^p{q,g} and ^^slgx- The former set is consistent with the integrability properties of the 
model, while the latter is more suitable for high-energy QCD as it reveals the Bose properties of 
the corresponding A^— gluon states. 



3. Baxter Q— operator 



Solving the Schrodinger equation (|2.8| ), we shall apply the powerful method of the Q— operator 
[pi][]. This method plays an important role in the theory of integrable models as it provides 
an alternative to the conventional Algebraic Bethe Ansatz. It is based on the existence of the 
operator Q that acts on the quantum space of the model and satisfies a finite-difference Baxter 
operator relations. In contrast with the ABA, the method of the Q— operator does not assume the 
existence of highest weight representation for the eigenstates and, as a consequence, it has a wider 
range of applicability. Both methods become equivalent if the eigenvalues of the Q— operator are 
restricted to be polynomials in a spectral parameter. In this case, the Baxter equations are 
reduced to the Bethe equations on the roots of these polynomials. It is not obvious, however, 
whether polynomial Q— operators furnish all relevant physical solutions to the Baxter equation, 
or equivalently the ABA method is complete. This turns out to be the case for the SL{2,M.) 
Heisenberg spin magnets ||21], |T^, while for the noncompact, SL{2,C) spin chain one has to go 
beyond the class of polynomial solutions. 

In this Section we shall construct a general (nonpolynomial) Q— operator for the homogeneous 
SL{2,C) spin chain that we shall denote as Q{u,u). This operator acts on the quantum space 
of the model V ... V, with V = V^'^'^\ and depends on two spectral parameters u and u. 
As we will see in a moment, for Q(u, u) to be a well-defined operator, these parameters have to 
satisfy the same condition ( p.21| ) as those for the i?— matrix. Following [^, we require that the 
operator Q{u, u) has to satisfy the relations 

• Commutativity: 

• Q — t relations: 

• Baxter equations: 



,u),Q{v,v)] = 0. (3.1) 

[tNiu), Q(m,m)] = [tjv(M), Qiu,u)] = 0. (3.2) 



tNiu)Qiu,u) = 


= {u + is)^(Q{u + i,u) + {u- 


-tsfQiu-i,u), 


(3.3) 


ijviu)Qiu,u) = 


= {u + is)^Q{u,u + i) + {u- 


-trsfQiu,u~t), 


(3.4) 



where t]\i{u) and tAr(ti) are the auxiliary transfer matrices defined in ( ^.47] ). 

According to ( |3.1|) , the Baxter Q-operator and the auxiliary transfer matrices, t]sf{u) and 
i]sf{u), share the common set of the eigenfunctions 

Q(M, U) ^p>,{g,g-}(£l, Z2, ..., Zn) = Q{q,q}{U, U) %,{g,g-}(£l, ^2 , •••, ^iV ) , (3.5) 
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which, at the same time, are the solutions to the Schrodinger equation (|2.8| ) The eigenvalues 
of the Q— operator, Q{u,u), satisfy the same Baxter equations, Eqs. ( |3.3D and (|3^) , with the 
auxiliary transfer matrices ( p.47| ) replaced by their corresponding eigenvalues. 



Our approach to constructing the Q— operator for the 5'L(2,C) spin chain is inspired by 
previous works |T2| , |I^, in which the Baxter Q— operator was constructed for the periodic Toda 
chain and the homogenous SL{2,M.) Heisenberg spin magnet. Namely, we shall represent the 
Q— operator by its integral kernel 



(3.6) 



and find the explicit form of the kernel by solving Eqs. ( p.lD -( p^ 



3.1. Baxter equations 



According to the definition ( |2.47| ), the auxiliary transfer matrix tjq{u) is a polynomial of degree 



A^ in the S'L(2,C) spin operators Sa^ . Using (|2.47|) , tN{u) can be expressed as the A^— order 
differential operator acting on the holomorphic coordinates Zk- Its substitution into the l.h.s. 
of ( p.3|) leads to a complicated holomorphic differential equation on the kernel Qu,u- One finds 
similar antiholomorphic equation from (|3.4| ). 

Instead of trying to solve the resulting differential equations we follow the approach developed 
in |l^, Q. It allows to find the exact solution to (|3.3| ) and ( |3.4D by using the invariance of the 
auxiliary transfer matrix t]s[{u) = trTjv('u) under local (gauge) rotations of the Lax operators 

Lk{u) ^ Lk{u) = M^n^iu) Mk+i , Tn{u) ^ fN{u) = M^^ Tn{u) M, , (3.7) 

where M^ are arbitrary 2x2 matrices, such that M^+i = Mi and det M^ 7^ 0. Let us choose 
the matrices M^ as 

with yi, ■■■, Hn being arbitrary gauge parameters. The matrix elements of the rotated Lax oper- 
ator ( ^.71 ) are given by 

[Lfc]ii =u + is + i{zk - yk+i)dzk , [Lk\22 = u - is - i{zk - yk)dzk , 

[Zfclai = K^k - yk)izk - yk+i)dz^ + (m + is){zk - yk) + {~u + is){zk - yk+i) , (3.9) 

while [Lfc]i2 remains unchanged. 

Let us now define the following function 

N 

y!^\z, y) = l[ [zk - yk]-'-'"- [zk - t/.+i]-^+^" , (3.10) 

fc=i 

where [z — yj-'^i'" = (^2; — y)~*^**'(^ — 1/)"'^=''*", ^at+i = yi and z = (i*i, ..., zjy). It depends on A^ 
auxiliary vectors y = {yi, ...,y^), as well as the spins (s, s) and the spectral parameters {u,u). 
The unique feature of the function F„ ^J* {z, y) is that it is annihilated by the off-diagonal matrix 
element of the gauge transformed Lax operator, 

[Lk],,Y^f\z,y) = Q (3.11) 



for k = 1, ..., N. Therefore, applying Y^^^'^ {z, y) to the matrix elements of the Lax operator L^, 
one finds that the latter takes the form of the upper triangular matrix. The product of such 
matrices can be easily calculated leading to the monodromy operator of the form 



u,u 



(«-^^)^^!:a 



(«+^^)^via 



(3.12) 



where in the last relation we used the explicit expression for the diagonal elements of the Lax 
operator ( |3.9| ). Taking the trace in the both sides of this relation, we obtain the following relation 
for the transfer matrix t^v = trT/v(M) 

t^{u; S^^\z)) Y^\z, y) = {u- zsf y!^,{z, y) + {u + zsf y!^,{z, y) . (3.13) 

Here, we indicated explicitly that the transfer matrix tN{u) is expressed in terms of the differential 
operators S^^\z) acting on the z— coordinates of the particles and defined in (|2.2| ). Remarkably 
enough, the relation (|3.13D takes the form of the Baxter equation (|3.3|) and, as we shall see later 



in this Section, it allows to construct the Q— operator. 

Repeating similar analysis for the transfer matrix in the antiholomorphic sector, iN{u), one 
can show that the same function (|3.10|) satisfies the antiholomorphic Baxter relation 



i^{u- S^'\z)) Yff\z, y) = {u- isf Fiji(z, y) + (n + isf ri;,5,(z, y) . (3.14) 

Notice that in order for the function Y^^;^^'{z,y) to be well-defined on the plane, the spectral 
parameters u and u have to satisfy the condition ( p. 211) . 



3.1.1. Properties of the F— function 

Before we proceed with constructing the Baxter Q— operator, let us discuss some properties of 
the function Yufi{z,y), which will be used below. 

Since the monodromy matrix Tf^iu) is related to the operator T]sr{u) by the gauge transfor- 
mation ( p.7|) , its matrix elements can be expressed, using (p.44|) , in terms of the operators An{u), 



..., Dn{u) and the gauge parameter yi. Substituting the resulting expressions into ( p.l2| ), we 
find that the function Yj'^^ (z, y) satisfies the following relations for arbitrary yk 

[ylB^iu) + yiiA^iu) - Dn{u)) - Cn{u)] Y!;f{z, y) = (3.15) 

and 

[A^{u) + y,B^{u)] Yjf {z, y) = {u - tsf FJ!f,(^, v) , 

[Dr,{u) - yiBr^iu)] Y^fiz, y) = {u + tsf Y^f^iz, y) . (3.16) 

We recall that the operators An{u), ..., D^iu) act on the holomorphic z— coordinates of A^— particles 
and do not depend on the gauge parameters yu- 
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Let us examine the relations ( |3.15|) and (|3.16|) in the hmit yi — > oo. At large yi, we use the 
definition ( p. 101) and expand the F— function as 



y!if{s,y) = [-yif {Ai%'^(zi, ...,iW|y2, -Jn) + 0(i/z/i, iM) 

with the leading term given by 

A (s,s) / -> -* I -• -* \ 



(3.17) 



^N-1 



[zi - 2/2J 



n [^fc - 2/d"'"'" i^k - vk+ir'^'"' 1 i^N - z/iv]-^-™ 



,fc=2 



(3.18) 

Then, one finds from ( 3.15 ) and ( 3.1(: ) that the function A^^^* , defined in this way, satisfies the 
following relations 

B,,iv)Ai:f{z\ y) = 0, A^iv)Ai:f = {u + ^sfA^^^, , D^iv^A^^ff = {u - ^sfA^fl^, . 

(3.19) 
Obviously, the same relations hold in the antiholomorphic sector. As we will show in Sect. 4.2, the 
function A^^^ (z | y) becomes a building block in the construction of the unitary transformation 
to the Separated Variables. 

3.1.2. The kernel of the Q— operator 

Taking into account the properties of the function F^/ , Eqs. ( p.l3|) and (|3.14| ), we look for the 
kernel of the Baxter Q— operator in the form 



Qu,ui^\w) 



.{s,s) 



d yYl^ '{z,y)ZR{y,w) 



(3.20) 



where (fy = (fyi...(Pyj\f and Zji is assumed to be a well-defined function on the plane, inde- 
pendent on the spectral parameters u and u. Substituting this ansatz into ( |3.(i| ) one finds that, 
thanks to Eqs. ( p.l3| ) and ( p.l4| ), the Baxter equations ( p.3| ) and ( p.4| ) are satisfied for arbitrary 
function ZR{y,w). 

To fix the form of the function Zfi{y, w) in ( |3.20| ) we require that the Q— operator has to 
commute with the auxiliary transfer matrices, Eq. ( p.2|) . It is convenient to represent the same 
condition in the form of the Baxter equations, Eqs. ( |3.3| ) and ( |3.4| ), with the transfer matrices and 
the Q— operator interchanged in the l.h.s. of these equations. Applying the operator Q{u, u) tj^{u) 
to an arbitrary test function and replacing the Q— operator by its integral representation, one 
integrates by parts to arrive at the relation analogous to (p.l8|) 



u,u)tN{u)^]iz) = / £w {tN{u;-S^'-'\w))Qu,ui^\w))'^{w) 



(3.21) 



where d'^w = (Pwi... (Pwn and the auxiliary transfer matrix in the r.h.s. is obtained from ( 2.47|) 
and (|2.44| ) by replacing the holomorphic SL{2, C) generators in the expression for the Lax oper- 
ators as Sa -^ —Sa * . In this way, the condition (|3.2|) on the Q— operator can be formulated 



as 



(-1)^ tM{-U; S^'-'\w)) Qu,u{z\ W) = {U + tsf Qu+i,u{^\ w) + {u - tsf Qu-^,u{^\ w) , (3.22) 
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where we took into account that the auxihary transfer matrix is a polynomial of degree N in 
the spectral parameter with the coefficient in front of u^ proportional to the (A^ — k)th power 
of the spin operators. Similar relation holds for the antiholomorphic auxiliary transfer matrix 



Let us look for the solution to (|3.22| ) in the form 

Qu,u{z\ w) = [a{s + iu,s- m)]^ Q^u-uiw \ z) 



(3.23) 



with the functions a(x) defined in (|A.3|) . Substituting this ansatz into (|3.22|) and changing the 
spectral parameter as m — > —u^ one finds that the function Qu,u satisfies the Baxter equation 

t^{u,S^^-'\w))QuAw\z) = {u + i{l-s))'' Qu+Lu{w\z) + {u-t{l-s)f Q^^i,^{iu\z) , (3.24) 



which coincides with (p.l3|) once one changes the spin as s -^ 1 — s. Therefore, the general 
solution ( p. 231) to the Baxter equation (|3.22|) and its antiholomorphic counterpart can be written 
as 

Qu,uiz\w) = [a{s + iu,s-iu)]^ / d'^y ZL{y,z)Y!:];;^''^~''\w,y) , (3.25) 

with Z^ly, z) being a well-defined function on the plane, independent on the spectral parameters 
u and u. 

Comparing two different representations for the kernel of the Q— operator, Eqs. ( p. 20 ) and 
( p.25|) , we notice that ( p.20| ) fixes the dependence of Qu,u{,z\ w) on the 2;— coordinates, whereas 
( p.25|) fixes its w— dependence. We find that ( p.20| ) and ( p.25|) become compatible provided that 
Zr is given by the following expressions 



TV 



N 



zUy.^) = \{ 



Wk-i - yk 



l2s~2 



Zp, (y, w) = Y\_ S'^{wk - yk) [wk - Wk+i] 



2s-l 



(3.26) 



fc=i 



k=l 



with wq = wn and wn+i = wi. Substitution of (|3.26|) into (|3.20|) yields two different expressions 
for the kernel of the Q— operator 



Qzmw) 



N 



d^y X{[zk - yfc]-^-™ [zk-i - 2/fc]-^+'" K-1 - Vkf'-^ 



(3.27) 



fc=i 



N 



and 

Qul{z\ilj) 



[a{2 -2s,s + %u, s - iu)7^f JJK - ^fc]'"'+'™ K " ^fc+i]'"'"™ l^k - Zk+i]'-'' 

k=l 



N 



n [zk - Wk]-'-''' [zk - wfc+i]-^+™ K - wk+i]''-' 



(3.28) 



k=l 



[a(l — 2s,s + iu, s — iu)/TT] 



N 



N 



d^ l[[wk - z/fe]^-'+*" K-i - ykY-'-''' [zk-i - yk]-'' 



k=l 



Here, the second relation in both equations is obtained from the star-triangle identity ( [A.8| ). The 
diagrammatical representation of the kernels Q'i^^l and Q^~^ is shown in Figs. || and |^, respectively. 
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1-s-iu 



^A(u,u) 




^ ^ 



^1 2s-1 ^ 2s-1 ^3 ^N 2s-1 ^1 



Figure 6: Diagrammatical representation for the kernel of the operator Q+(n,u) defined in 
Eq. (13:271) ■ v4(m, m) = [vr a(2 - 2s, s + m, s - iu)f . 



1 1-2s ™2 1-2s ™3 ™N 1-2s ™1 




^ B(U,U) 




Figure 7: Diagrammatical representations for the kernel of the operator Q_{u,u) defined in 



Eq. (lOSi) . B{u, u) = [a{l -2s,s + iu, s - iu)/ 



71 



iN 



Using this representation, one checks that the expressions ( |3.27| ) and ( ^.28| ) match into Eqs. ( |3.2CI| ) 
and ( p.25| ) simultaneously. 

The obtained expressions for the kernels, Eqs. ( p.27| ) and (|3.28|) , define two different operators, 
Q+(u, u) and Q_(it, u). By the construction, they commute with the auxiliary transfer matrices, 
Eq. (|3.2|) , and satisfy the Baxter equation (|3.3|) and (|3.4| ). For the kernels of these operators, 
Eqs. ( p. 27] ) and ( p.28| ), to be well-defined on the plane one has to require that, similar to the 
i?— matrix, the spectral parameters u and u have to satisfy the condition ( p.21|) . Finally, in order 
to identify Q+{u,u) and Q_(u, u) as the Baxter Q— operators we have to show that they fulfil 
the commutativity condition (|3rT|) . 



3.2. Commutativity condition 

Let us show that the operators Q+{u,u) and Q-{u,u) with the kernels defined in ( |3.27| ) and 
( p.28|) commute with each other for different values of the spectral parameters 



and satisfy the following exchange relations 



u,uu%P-[v,v 



S^+l") "^ 



(u, u) Q^v, v) = Q+{v, v) Q_(u, u) 
{u, u) Q+{v, v) = Q_{v, v) Q+(u, u) 



(3.29) 
(3.30) 



We start with the product of two operators 

[Q_{v,v)Q+{u,u)]{z\w) 



d yQv,l{z\y)QlUy\w) 



(3.31) 



and substitute them by the corresponding Feynman diagrams. It becomes convenient to use the 
left diagram in Fig. ^for the operator Q_(f,'u) and the right diagram in Fig. ^for the operator 
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Figure 8: Diagrammatical representation of the kernel Xu,u;v,v{z,w) defined in Eq. ( p.32| ) as 
a periodic chain of rhombuses. The rightmost vertex is identified with the leftmost one. To 
prove the commutativity relation ( |3.35| ), one inserts two lines with opposite indices into the left 
rhombus between the points zi and wi. 



Q+(m, m). Gluing together two sets of the triangles we notice that their common lines with the 
indices (1 — 2s) and (2s — 1) annihilate each other and the resulting Feynman diagram takes the 
form shown in Fig. ^. The corresponding Feynman integral can be written in a simple form by 
introducing notation for the following function 



Xv,v;u,u{^\ w) = [a{s + iu,s- iu)] / d yY^y'{y, z) Y\X iv, ^) 



'(«.«)/ 



-(1-5,1-5), 



with the y— functions defined in ( p.lOl ). In this way, one arrives at 



N 



P_(t;, v) Q+(u, u)] {z\ w) = [7ra(2 - 2s)] X^^jj.u,uiz\ w) 



(3.32) 



(3.33) 



Following the same steps, we now calculate the product of the same operators but in an opposite 
order, [Q^{u,u)Q-{v,v)]{z\w). This time we replace two operators by the left diagram in 
Fig. P and the right diagram in Fig. 0, respectively, and obtain two sets of the star-diagrams 
glued together through common vertices. Integration over the position of these vertices can be 
easily performed using (|A.7|) and it gives rise to the 5— function connecting the centers of the 
star- diagrams. As a result, one arrives at the diagram, which is similar to the one shown in 
Fig. P, leading to 



N 



[Q4u,u)Q4v,v)]{z\w) = [7ra(2 - 2s)]'' Xu,u;vAz\w) 



(3.34) 



Comparing ( p.33|) and (|3.34|) we notice that their r.h.s. differ from each other by interchanging 
the spectral parameters. Then, the commutativity of the Q+(m, tt) and Q^{v,v) as well as the 
relations ( |3.3CI| ) follow from the following symmetry property of the X— function 



-^u,u;v,v\Z\'W ) -^v,v ■,u,u\^\'W ) ■ 



(3.35) 



The proof of ( ^.35| ) is based on the uniqueness relations, Eq. ( |A.8| ), and it can be carried out 
using the diagrammatical representation of the function Xu^u;v,v- It also relies on the permutation 
identity shown in Fig. |^. To verify this identity it is sufficient to turn the "unique" triangles in 
the both sides of the relation into unique stars and check that the resulting diagrams coincide. 
Turning to (|3.35|) , we use the diagrammatical representation of the function Xu^u;v,v and insert 
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a(s+iv) a(s- iv) 



s-iu 




1-s+iv 



s+iv 



1-s-iu 




3 
I 
> 



a(s+iu) a(s-iu) 



Figure 9: The permutation identity. 



two additional propagators with the indices i{u — v) and —i{y — u), respectively, into one of the 
rhombuses as shown in Fig. H. Since the sum of the indices vanishes, this transformation does 
not change the function. Subsequently applying the permutation identity (see Fig. |^), we move 
the line with the index i{v — u) to the right of the diagram until it returns to its initial position 
and annihilates the line with an opposite index, thanks to periodic boundary conditions along 
the chain of rhombuses in Fig. ^ In this way, one arrives at the initial diagram, in which the 
spectral parameters u and v are interchanged, thus proving the relation (|3.35| ). 

Let us now turn to the first two relations in ( |3.29|) and examine the product Q+(m, u) Q+(f , v). 
Representing the first operator by the right diagram in Fig. ^ and the second one by the left 
diagram, one obtains the following Feynman integral 



N [ ^ 

^ k=i 



Zk+lY ^^Wk 



Vk+i 



\2{l-s) 



(3.36) 



Due to the symmetry property ( p.35|) , the r.h.s. of this relation is symmetric under permutation 
of the spectral paramaters and, as a consequence, the operators Q+(m, u) and Q+(w, v) commute. 
The proof of the commutativity of the operators Q-(m, u) and Q-(v, v) goes along the same lines. 
Finally, Q+(m,'u) and Q_(m,m) are the 5'L(2,C) invariant operators 



[Q+(m,m),S'q 



{u,u),Sc 







(3.37) 



^^^ 



with Sa = ^fc=i Sa being the total SL{2, C) spin of the system. The same property implies 



that the kernels of the Q— operators have to transform under the SL{2, C) transformations (2~3) 



as 



N 



Wt::;' 






Z W 



Ulcwk + df-^'[czk + df' g2(z I w) . 



(3.38) 



^fc=i 



This relation can be verified using the explicit form of the kernels, Eqs. ( |3.27| ) and ( ^3.28] ). 



3.3. Properties of the Q— operator 

In the previous Section we have constructed two different Baxter Q— operators, Q+(u,u) and 
Q_(ii, u). Each of them satisfies the defining relations, Eqs. ( |3.1j ) - ( p.4| ), as well as the additional 
relations ( p.30|) and (|3.37|) . Let us show that these two operators are conjugated to each other 
with respect to the scalar product (|2.9|). 
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3.3.1. Conjugated Q— operator 

The kernel of the operator Q^{u,u) conjugated to the Baxter Q— operator is defined as follows 

(ti,u)]^^] {zi,...,zn) = / (fw {Q^^a{wi,...,WN \zi,...,ZN)y '^{wi,...,wn). (3.39) 



Calculating [Q+{u,u)]^ , one substitutes its kernel, Eq. (|3.27|) , into the r.h.s. of ( p.39| ) and notices, 
using ( p.l2| ), that ((5i^^(w | z)) turns out to be proportional to the kernel Q^^'^. (^| w) defined 
in ( p.28| ). This leads to the following relation between the corresponding Q— operators 

[Q+(m, u)]^ = (a(l - 2s, s + m*, s - iu*)/ny^ Q-iu*, u*) . (3.40) 

The inverse relation looks like 

[Q_(u, u)]^ = [(a(2 -2s,s + iu*, s - iu*) tt)"^ Q+{u*, u*) . (3.41) 

Since the operators Ql.{u,u) and Q^(u*,u*) differ from each other by c-valued coefficient func- 



tion, they commute with each other for different values of the spectral parameters due to ( |3.29|) 



[ [Q±{u, u)]\ Q^iv, v)] = [ [Q±(m, u)]\ Q±{v, v)]=0. (3.42) 

Additional properties of the operators Qj.{u,u) follow from conjugation of the corresponding 
relations for the operators Q±{u,u). In particular, the conjugated Baxter equations ( p.3|) and 
look as 



tNiu) [Q±(r , M*)]t = {u- 2(1 - s)f[Q^iu* - t, u*)]^ + iu + 2(1 - s)f[Q±iu* + 2, u*)]^ , 

tNiu) [Q±(r, M*)]t = {u- 2(1 - rs)f[Q±iu*, u* - 2)]t + {u + 2(1 - rs)f[Q±iu\ u* + 2)]+ . 

(3.43) 

Here, the spectral parameters u and u are arbitrary complex numbers satisfying the condition 



The relations ( |3.4CI|) - ( |3.43|) can be further simplified if one imposes the additional condition 



on the spectral parameters, u = u*, leading to ( ^.39| ). As we will show in Sect. 4.4, it is for these 
values of the spectral parameters that the Baxter operator Q{u, u*) enters into the expression 
for the eigenstates of the Hamiltonian in the representation of the Separated Variables. Notice, 
however, that in the r.h.s. of the Baxter equations ( p.43| ) the holomorphic and antiholomorphic 
arguments of the Q— operators are shifted by ±2 independently and, therefore, one can not obtain 
from ( 3.43 ) a closed system of equations on the operator Q(m, u*). 



3.3.2. Quasimomentum 

The Baxter Q— operator commutes with the transfer matrices of the model and, therefore, it 
could serve as a generating function of the Hamiltonian and the integrals of motion. It is well- 
known from the theory of the compact spin magnets that the main observables (like Hamiltonian, 
quasimomentum, etc) are determined by the behaviour of the Baxter Q— function at the special 
values of the spectral parameter u = ±is [^. As we will show below, the same formulae hold 



for the noncompact spin magnet. 
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Let us examine the operator Q+ {u, u) with the spectral parameters given by u = ±zs and 
u = ±is. Using the diagrammatical representation of this operator, Fig. |^, we find that the lines 
with the indices s ±iu disappear and the integrals over the centers of the star diagrams take the 
form of (|A.7|) and give rise to the 5— function 



N 



N 



Qa(i'i^)=p;^n^^'^( 



Zk - Wk) , 



Q 



(+) 



.Z\W] 



-N 



l[6^^\zk+i-Wk), (3.44) 



fc=i 



fc=i 



where ps = l/(7r^a(2s, 2 — 2s)) = (n^ + 4z/^)/7r^. Going over to the operators one gets 

Q+iis,is) = p;^ X 1, Q+(-^s, -is) = p;^ x P, (3.45) 



where P is the operator of the cyclic permutations defined in (|2.58|) . Putting u = ±is and u = ±is 
in ( |3.40| ) and taking into account ( |3.45|) , one finds that the a— factor in the r.h.s. of ( |3.4CI| ) vanishes 
and, as a consequence, the operator Q_{±is*, ±is*) is divergent. A careful analysis shows that 



.{-is* + e, -is* + e) = 1 X 



Uvr 



\N 



ns — e,is 



P-^ X 



[m 






■[^ + 0{e)] 



cN 



■[^ + 0{e)] 



(3.46) 



as e — > 0, so that the Q_— operator exhibits the iV— th order pole. The detailed analysis of the 
analytical properties of the operators Q±(m,'u) will be performed in Sect. 3.5. 

Using (|3.45|) and (|3.46|) , one obtains the following relations for the operator of quasimomentum 
0{q,q) defined in Eq. (|2.65|) 



e 



-z In P = i In ■ 



Us,«s 



iln ■ 



-IS, —IS] 



-is* + e, —is* + e) 



e-^O 



(3.47) 



We recall that the eigenvalues of this operator have the general form ( 2.66 ) and satisfy the relation 

3.3.3. Properties of the eigenvalues 

As we have seen in Sect. 2.4, the Hamiltonian of the model is invariant under the cyclic and 
mirror permutations generated by the operators P and M, respectively. Let us now examine the 
action of these operators on the Q— operators. Since the Q— operators commute with each other 
for different values of the spectral parameters, it follows from ( p.45|) that 



i±[u,u) 



0. 



(3.48) 



The same property can be easily understood from Figs. ^ and |^ - the diagrams representing the 
Q— operators remain unchanged if one performs simultaneously the cyclic permutations of the 
z— and w— vertices. 

The tranformation properties of the Q— operators under mirror permutations become more 
complicated and can be expressed as 



»+(-M, -n) =MP"^ Q+(m,m) 



.{-u,-u) 



:p 



[U,U 



(3.49) 
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To see this, one applies the operators M and P to the diagrams shown in Figs. ^ and |^ and finds 



that the diagrams describing the r.h.s. of the relations (|3.49| ) look like mirror images of those 
shown in Figs. ^ and |^. To bring them back to the original form, one fiips the whole diagram along 
the vertical axis and changes simultaneously the sign of the spectral parameters as indicated in 
the l.h.s. of ( |3.49| ). One can verify ( p.49| ) by putting u = —is and u = is* in the first and the 
second relations, respectively, and substituting the Q— operators by their expressions ( p.45| ) and 

Applying the both sides of (|3.49|) to the eigenstate '^p,{g,q} one finds the relation between the 
corresponding eigenvalues of the Q— operators defined in (|3.5|) 



-M, -m)|^ 



p,{g,g 



}) 



:p" 



[u,u 



%,{-g,-g}) 

^{-l-.-}(^'^)I^P^-{M}), (3.50) 

q) = {{—l)^qk, {—l)''qk)- Here, in the second relation we used Eq. ( |2.64D and in 



-ie{-q~q) ^(+) 



where (— g 

third one Eqs. (|2.65| ) and ( |2.67|) . Then, taking into account (|2.65|) we find 



Q 






-u, —u) 



MM) 



Q^+^ -Au,u) 
^{-q-q}^ ' ' 



In a similar manner, we find from the second relation in (|3.49|) 



^U}(""'"^) 



-ie{q,q) n(~) 



QV' 



{-9,-g-}V", 



u,u] 



Let us now consider the ratio of the operators 



p:""' X 



(3.51) 



(3.52) 



(3.53) 



with ps = {nl + 4z/^)/7r^. According to ( p.30| ), the operator W defined in this way does not 
depends on the spectral parameters. The normalization factor is chosen in ( |3.53| ) in such a way 
that the operator W is unitary 

U{u*,u*)yQ+{u,u) 



W^W 



P. 



N 



{u*,u*)\ Q_{u,u) 



1. 



^ 



p,{q,q} 



Qi'Wq,q ^ 



p,{q,q} ■ 



(3.54) 



Here, we used the fact that the operator W does not depend on the spectral parameters and 
applied the identity (|3.40|) . Since the eigenstates \l/p-^{g^g}(2') diagonalize the operators Q±, they 
also diagonalize W. Then, we obtain from ( |3.53| ) and ( p.40| ) 

^{M}("'^) _.....u.„ , ... . ^u)e^-^f, (3.55) 



Q 



(+) 

{q,q} 



u\ u 



e™"-' \a{s + iu, s 



with e*"^' = [a(l - 2s, 2 - 2s)Y'^, or ^, = arg[r(l - n, + 2zi^,)r(-n, + 2ii^s 



We shall use the relation (|3.55| ) below to fix the phase of Q^ -y{u,u*) at u = u* and to 
calculate the residues of Q\VqAu,u) at its poles. For instance, putting u = ±(is* — e) and 

u = ±{is* — e) in ( |3.55| ) and applying ( p.45| ), we calculate the residue of Q 
pole in e as 



r -| at the A^— th order 



with Cs = iTips 



QHq}(-'^* + ^^-'^' 


' + 


A — JWq,q S 


1 + 


0{e)] 


Q[:'d^-^* 


— e,is* — 


e) = 


, Q-iO,. 


■<?!:•„( 


-is* 


+ e,- 


'^' = ^[in^s 


+ 4z/2)/7r4] 


-1/^ 


' and . 


3 = 1- 


S*. 





(3.56) 
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3.3.4. Reduction formulae 

The operator Q+(m, u) acts on the quantum space of the system Vi ® ... ® Vn, with Vk = V^^'^\ 
and depends exphcitly on the number of particles A^. There exist the relations connecting the 
Baxter operators defined for the different number of particles. 

To obtain the reduction relations, one examines the kernel of the operator, Q^uM^"^)' ^^ ^^^ 
limit when two particles with coordinates, say, zn and zi approach each other on the plane, 
\zn — Zi\ — > 0. Using the diagrammatical representation of the kernel - the left diagram in Fig. |^, 
one finds that in this limit two lines connecting the center of the left unique star diagram with 
the points z^ and zi merge and produce a single line with the index 2s. Then, the resulting 
integral over the center of the star takes the form ( |A.9| ) and gives rise to S{z]\f — wn) 

QuU^i,--,zn-i,zn\wi,...,wn-i,wn) (3.57) 

= 7r^a(2s,2(l - s))6{zn - wn) Qu,Uzi, ..., zn-i\wi, ...,wn-i) ■ 

Here, the Q+— operator in the r.h.s. acts on the quantum space of A^ — 1 particles, Vi ... ® Vat.i. 
One can continue the reduction procedure by putting any pair of the nearest neighbors atop of 
each other, Zk = Zk+i- Each time the number of particles is reduced by 1 until one reaches the 
limit zi = Z2 = ■■■ = zm when the iV— particle Baxter Q— operator is degenerated into the identity 
operator ~ Ilfc '^(-^i ~ "^fc)- Similar reduction formula holds for the kernel of the Q_— operator. 



Eq. (|3.28|) , but in this case one merges the right arguments of the kernel 



Q\,}i{zi, ...,zn-i,zn\wi, ...,wn-i,wn) ^ ^ (3.58) 

= 7ra(2(l ~s),s + iu,s- iu) 5{zm - wn) <5^(^i, •••, zn-i\wi, ...,wn-i) ■ 

Here, the additional factor comes from the 5— factor in the r.h.s. of Fig. |^. 

Let us now put zi = zjy in the both sides of the eigenproblem ( |3.5|) and apply the relation 

(iin 



7r2 



QSi\u' m) ^p,{g,g-}(i'l, ..., ZN-1, ^l) = ^2 + 4^2 



X d wQl^l,{Zi,...,ZN-l\Wl,-..,WN-l)'^p,{q,q}{Wi,...,WN^l,Zi),{3.59) 



where the superscript (A^) in the l.h.s. indicates that QqJ{u, u) is the eigenvalue of the A^— particle 
Baxter operator. Using the property of completeness of the solutions to the eigenproblem (^]5|), 
we may expand '^^^{q^q^iwi, ...,wn-i,zi) over the eigenstates of the (A^ — 1)— particle Baxter 
operator, ^!^|„Vg/|(wi, ...,wn-i) with the expansion coefficients depending on zi. After its sub- 
stitution into (|3.59|) , the integral in the r.h.s. is replaced by the sum over the eigensvalues Q^^Zp . 
Then, comparing the z— dependence of the both sides of the resulting relation one arrives at 
overcompleted system of equations on the expansion coefficients. Its consistency requirements 
provide the quantization conditions on the integrals of motion g^ and g^ [0. Solving the sys- 
tem one finds the expansion coefficients and, as a byproduct, reconstructs the expansion of the 
eigenfunctions and the eigenvalues of the A^— particle Baxter equation, Qq^^{u,u), over those 
corresponding to the system of (A^ — 1)— particles [p^] . 
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3.3.5. Eigenvalues at iV = 2 

The eigenproblem for the Baxter operator Q+(m, u) can be solved exactly for the system of A^ = 2 
particles. In this case, the eigenfunctions of Q+(u, u) are given by ( |2.55| ). They are parameterized 
by the pair of the SL{2, C) spins {h, h) and by the total momentum p, or equivalently the center- 
of-mass coordinate zq. Their substitution into (p.5|) yields 



(fwi(fw2 Qufii^l, Z2\W1, W2) ^fc,h(u^l - Zq, W2 - Zq) = Qh,h{u, u) ^/i,h(^l - ^0, ^2 - Zq) , (3.60) 

where the function '^f^j^{zi,Z2) was defined in ( 1^.551 ). Following the diagrammatical approach, 
we represent the integral entering the l.h.s. of this relation as the Feynman diagram shown in 



Fig. |TO|a. To arrive at this diagram we replaced the kernel of the Q— operator by the right diagram 
in Fig. H at A^ = 2 and represented the eigenfunction '^f^f^{wl — Zq, W2 — ^o) as a (nonunique) 
triangle with the vertices at the points Wi, td2 and zq. 

The both sides of the relation ( |3.60| ) depend on two vectors Zi — zq and 22 — Zq that we may 
choose to our convenience. One possible choice could be ii — > 00, io = and Z2 = Z2 = I. In this 
limit, four lines connecting the point Zi with the rest of the diagram can be effectively removed 
and the resulting diagram takes the form shown in Fig. p!0|b. The corresponding Feynman integral 
defines the eigenvalues of the Baxter operator Q+{u, u) at N = 2 



(tt a(2 — 2s,s + iu, s — iu)) 



2 



U) T '[u,u) = I a Wid W2 7 — r; 7—. — r^ —ft ttttt rrrrT ttt^ — r • (o.ol) 



One can find another (equivalent) representation for the same eigenvalue by choosing in ( |3.6U| ) 



io — > 00, 2*2 = and Zi = Zi = 1. In this case, one can remove two lines attached to the point 



Zq. The resulting diagram takes the form shown in Fig. |10p. By the construction, Q 7 {u,u 



^(7V=2), 



satisfies the Baxter equations, Eqs. (|3.3| ) and (|3.4| ) for N = 2, with the auxiliary transfer matrices 
given by ( p.47|) 

t2iu) = 2u'^ ~ h{h - 1) + 2s{s ~ 1) (3.62) 

and i2{u) is given by a similar expression in the antiholomorphic sector. 

The two-dimensional integral in (|3.61| ) has a striking resemblance to expressions for the cor- 
relation functions in two-dimensional conformal field theories and, therefore, can be evaluated 
using the powerful technique developed in p5[. The resulting expression for ( p.61| ) is given by 



the sum over the product of holomorphic and antiholomorphic "conformal blocks" 
^hl"^^(^'^) = c{s,h) 



,{Ar=2)/ ^_,^ „^^ ^^T{1- s - iu)T{l - s + iu) 



r(s — iu)T{s + iu) 
xlQoiu) {Qoi-u*)r + (-1)"" Qoi-n) {Qo{u*)r \ , (3-63) 



where h = {1 + nh)/2 + ivh, h = 1 — h* and Qq{u) is given by one-dimensional contour integral 
that can be expressed in terms of hypergeometric series as 

nr\ f2s-h,2s-l + h,s + m\ 

Qoiu)=sF2i ^^ 1 . (3.64) 



The constant c(s, h) in ( |3.63| ) does not depend on the spectral parameter and it is fixed by the 



normalization condition (|3.45| ), Qi^j^{is,is) = vr^/(n^ + 4//^ )^. 
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(a) 



Figure 10: The eigenvalue of the Baxter 1 
representation; (b) reduction ai Zq = 1, Zi 
and Z2 = 0. 



(b) 



(c) 



-operator at A^ = 2: (a) general diagrammatical 

• oo and Z2 = 0; (b) reduction at zq = oo, zi = 1 



As a matter of fact, the expression ( ^.641 ) has already appeared in the analysis of the Baxter 
equation for the noncompact spin chains. It was introduced for the first time in [^ ^ as 
the solution to the holomorphic Baxter equation for the SL{2, C) magnet of spin s = 1 and 
s = 0. Later, the same expression was identified |1^ as the solution to Baxter equation for the 
Heisenberg SL{2,'R) magnet of arbitrary spin s and the total spin of the system, h — 2s, being 
nonnegative integer. In the latter case, Qo{u) is reduced to a polynomial of degree h — 2s in u, 
which turns out to be identical to the Hahn orthogonal polynomial. 

The detailed analysis of the properties of the obtained expressions, Eqs. ( p.63| ) and ( p.64| ). 



and their generalization to higher A^ will be presented elsewhere p4 . 



3.4. Relation to the Hamiltonian 



There exists the relation between the Hamiltonian of the model, ( |2.42| ), and the Baxter Q— operator. 
It allows to replace the original Schrodinger equation ( |2.8| ) by the spectral problem (|3.5| ) and ex- 
press the energy Ej^j in terms of the eigenvalues of the operator Q+(m, u). 

The above relation follows, in its turn, from a more general relation between the 5'L(2,C) 
transfer matrices of arbitrary spin, T^'^°'^°''{u,u), defined in ( p.40|) , and the Baxter Q— operator 
(see Appendix B for details). The simplest way to establish this relation is to compare the 
Feynman diagrams describing the transfer matrix and the product of Q— operators. Figs. ^ and 
1^, respectively. It is easy to see that, up to redefinition of the spins and the spectral parameters, 
two diagrams are identical. In this way, we obtain the following relation 



T 



N ^ 



X{u + i{l — so),u + i{l — so);u + iso, u + iso) 



(3.65) 



which is valid for arbitrary 5*17(2, C) spins (sq, sq) and the spectral parameters u and u. Here, 
the notation was introduced for the operator X(ti, u; v, v) with the kernel Xu^u;v,v{z, w) defined in 
( ^.32| ). Taking into account the relation between the X— operator and the Baxter Q— operators, 
Eqs. (|3.33| ) and (|3.34| ), we obtain the expression for the transfer matrix in terms of the operators 



T 



(so,so)( 
N ^ 



p'l'"'''\u,u 



{u,u) X Q_{u + i{l - so),u + i{l - so))Q+{u + iso,u + iso) 
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= p^^''''\u, u) X Q_(m + ISO, u + iso) Q+(m + i{l - So), u + ^(1 - So)) ,(3.66) 

where the spectral parameters satisfy the relation ( |2.21| ) and the normalization factor p^^^'^^^u, u) 
can be calculated from Eqs. (|2.25| ) and ( p.32|) as 



p. 



(so, so) I 



U, U) 



1 


a{so 


- s + iu, s — So - 


- iu) 


_7r3a(2 - 


■ 2s) a{s - 


- so + iu,so + s - 


- iu) 



N 



(3.67) 



The following comments are in order. 

According to (|2.42|) , the Hamiltonian of the model is given by the logarithmic derivative of 
the fundamental transfer matrix T*-*'*-'(u, u) at m = u = 0. Since the spectral parameters have 
to satisfy the condition (|2.21|) , to calculate the Hamiltonian it becomes sufficient to analyze the 
transfer matrix for real values of the spectral parameters u = u = u*. In this case, one puts 
So = s in ( p.66| ) and substitutes the Q_— operator by its expression in terms of (Q+)^— operator, 
( p:iDO , to find (see also ([B^)) 



T%^'\u,u) 



PqW I 

.2N 



(u — is,u — is)] ' 



(u + is,u + is) 



(3.68) 



PqU~" [(ii+('U — is* , u — is*)] Q+('U + is* , u + is*) 
with s* = 1 — s and the normalization factors 



PoW 



a{l — 2s + m, 1 — iu) 
a(2 — 2s — iu,l + iu) 



N 



Pc 



Pc 



nl + Avl 



2\ N 



TC^ 



One verifies, using (^), that T("'"")(0,0) = (-l)^p-^ Together with (|2:i2| ) this leads to 
(— l)^exp (— i^jv ^ iuTCiy + C^(^^))- We find the expression for the Hamiltonian 



N 



U, U] 



r 

Hn by substituting (|3.68| ) into this relation 
d 



au 



[u + ts,u + is] 



u=0 



'■ d , 
I—- m( 

du 



(u — is,u — IS] 



M=0 



where the notation was introduced for the additive normalization constant 



^N = i— In Pq{u) 



u=0 



2NRe [ip{2s) + ^(2 - 2s) - 2^(1)] . 



(3.69) 



(3.70) 



Eq. ( |3.69| ) establishes the relation between the Hamiltonian and the Baxter Q— operator. We 
recall that the latter satisfies the Baxter equations, Eqs. (|3.3| ) and (|3^) , and depends on the 
integrals of motion {q,q}. Then, Eq. ( ^.69] ) provides the explicit form of the dependence ( p^.54| ). 
Applying the eigenstate "^piq^q} to the both sides of (|3.69|) we calculate the energy as 



EN{q, q) = Sn + i{\nQg^g{is, is))' + i{\vL [Qq^q{- 



-IS, —IS 



)]*)': 



(3.71) 



(+) 



u,u] 



where prime denotes derivative with respect to the spectral parameter and Qq^q{u, u) = Qt^U 
stands for the eigenvalue of the operator Q+ («,«). We can further simplify this expression by 
using the properties of the eigenvalues ( p. 511 ) 



EN{q, q) =eN + «(ln Qq^q{is, is))' - i{\n [Q„q^_g{is, is)]*)' ^ 
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(3.72) 



where {—q,~q) = {{—l)'^qk, (— l)'^gfc)- Taking into account that EN^q.q) is invariant under the 
replacement (g, q) -^ (— g, — g) (see Sect. 2.4) one gets from the last relation the expression for 
the energy 



EN{q,q) =£n 



Im — In [Qq,q{u + is,u + is) Q-q-q{u + is,u + is)] 



u=0 



(3.73) 



which is explicitly real. This expression was obtained from the first relation in ( |3.68| ). Starting 
from the second relation in (|3.68| ), one can obtain another representation for the energy 



ENiq, g) = - Im — In [u^^Qq^q{u + is*, u + is*)Q^q^^q{u + is*,u + is' 



du 



u=0 



(3.74) 



We notice that, according to ( |3.56D , the Q— function entering this expression has the A^— th order 
pole at M = 0. The factor m^^ compensates this pole and, as a consequence, the energy is 
determined by the subleading 0{e) terms in the r.h.s. of (|3.56|) . 

One can verify |2^ that at A^ = 2 the substitution of ( p.63| ) into Eqs. ( p.73| ) and ( p.74| ) leads 
to the well-known expression for the energy (p.56 ). 



3.5. Analytical properties 

Let us examine the analytical properties of the Baxter operator, Q+(u,u), and its eigenvalues, 
Qq^q{u,u), as functions of the spectral parameters u and v,. We remind that the spectral param- 
eters have to satisfy the condition i{u — u) = n with n being integer. It becomes convenient to 
parameterize their possible values as m = —in/2 + A and u = in/2 + A with A being arbitrary 
complex. The operator Q+(A — in/2, A + in/2) becomes a function of integer n and continuous 
complex A. 

In this Section we shall determine analytical properties of Q+(A — m/2, A + in/2) on the 
complex A— plane for fixed n. To start with, we consider Eq. ( |3.5| ) and substitute the Q— operator 
by its integral representation with the kernel given by the first relation in (|3.271 ) 



N 



Vk-l] " "'" [Zk-l - Vk-l] '+™%,{5,5}(y) = Q{q^q}{u,u) ■^j;,{q,q}{z) 



(3.75) 



fc=l 



where the notation was introduced for the function 

N 



*P,{g,5}(^) 



I Swm 



ni 



P,{M]\w) I \[Wk-yk 

k=l 



\2s-2 



nt + Aui 



TT^ 



-N/2 



[m,^,{q£ (y) . (3.76) 



This relation is well-known as the SL{2, C) intertwining transformation with U being the corre- 

1T7| . It maps the state 'i!^^''^\z) belonging to the space V" ... V^, 



■s,l-s)i 



Z . 



sponding unitary operator 

with V = V^^''^\ into the unitary equivalent state \1/^^" 

The singularities of the eigenvalues Q{g,g}(A — in/2, A + in/2) in A are in one-to-one corre- 
spondence with divergences of the integral in the l.h.s. of ( p.75| ). The latter originate from two 
different integration regions: \yk — Zk+i\ -^ and \yk — Zk\ — > 0. In the first region, one expands 
the integrand in ( p.75[ ) in powers of Wk = yk — Zk+i, assuming that \l/p'^{g^g}(y) — \l/p',{g,g}(z) 



mu — l — mi. — 1 
UJi, 



and calculates the l.h.s. of ( p.75|) as the product of integrals of the form 



(f'Wk wl 



-'""wl 



-» m.- 



1 — mi. — 1 



27i6 



\w,,\<e 



s+iu—mk, s+iu—rrn; 

s + iu — rrik 



(3.77) 
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Figure 11: The position of poles (fat points) of the Baxter operator Q+(A — in/2, A + in/2) in 
the complex A— plane at n = 0. The SL{2, C) spin of particles, s = (1 + ns)/2 + ii/g, is chosen 
in such a way that n^ = and i/^ > 0. Two squares correspond to the values of the spectral 
parameters {u,u) = {±is,±is). 



with mfc and fhk being positive integer. Thus, integration over each ijk ~ Wfc+i in ( p.75| ) brings 
a simple pole in the spectral parameter located at s + iu — m^ = s + iu — fhk = 0. Taking 
their product we find that the resulting A^— fold integral over yk in the l.h.s. of ( 3.75| ) acquires 
the poles in the spectral parameter A of the order not higher than A^. Similarly, calculating the 
contribution from the second region y^ ~ Zk we find another (infinite) set of poles located at 
s — iu — ruk = s — iu — fhk = 0. 

Combining together the contribution from both regions, we conclude that the eigenvalue of 
the Baxter Q— operator, Q{q^qy{u,u), is an analytical function on the complex plane except the 
infinite set of points 



{ 



u^ 



lis 



m , u; 



i{s-m)} ; { 



u,. 



-i[s 



m] 



u.. 



-i[s 



m)} 



(3.78) 



with m, ffi = 1, 2, ..., at which it has the poles of the order not higher than A^. Using the above 
parameterization of the spectral parameters, u = —in/2 + A and m = in/2 + A, we obtain that 
the poles can be enumerated by a pair of positive integers (m, m) 



n; 



Tyfis + TTi 



m 



) ) K 



T\i^s + i- 



,m + m — 1 



(3.79) 



with s = (1 + ns)/2 + ii^g and ?ti, m = 1, 2, ... . Then, for arbitrary integer n = rim,™, the function 
Q{q^q}{X — in/2,X + in/2) has the iV— th order poles situated along the imaginary axis in the 
A— plane at the points A = Xm^m as shown in Fig. |ll|. These two infinite sets of poles lie on the 
different sides from the real axis and do not overlap. Since | ImA^^| > 1/2, the eigenvalues of 
the Baxter Q— operator do not have poles within the strip —1/2 < IniA < 1/2. 

As was shown in the previous section, different observables of the model - the Hamiltonian 
and quasimomentum, Eqs. ( p.73| ), (|3.74| ) and ( p. 47] ), are expressed in terms of the Baxter operator 
Q+ and its derivatives defined at the special points {u = ±is, u = ±is) and {u = ±is*, u = ±is*). 
In the former case, the points belong to the region of analyticity of the Q+— operator (see Fig. |Tl|) 
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and, therefore, the above operators are well-defined. In the latter case, the A^— th order pole of 
the Baxter Q+— operator is compensated by the factor u^^ entering the second relation in (|3.68|) . 

Analyzing the singularities of the Q+— operator we have used the integral representation of 
the kernel defined by the left diagram in Fig. p. One could use instead another representation, 
corresponding to the right diagram on the same figure. In this way, one finds that the Feynman 
integrals have poles aXu = ■iii{s + m) and u = ±i{s + m) with m and m positive integer, that are 
different from ( p.78| ). In addition, one has to take into account that the prefactor A{u,u) also 
generates singularities. One finds that A{u, u) produces the poles at the points (|3.78| ) and, at the 
same time, it vanishes at m = ±i(s + m) and u = ±i{s + rh) , thus compensating the singularities 
of the Feynman integral. 

Summarizing the analytical properties of the Baxter operator, we expand Q+(m,'u) into the 



sum over poles (|3.79| ) 



_(A-m/2,A + m/2) = "^ S. 



n."m,™ "*'™ 



m,m=l 



1 



+ 



zEj 



+ 6, 



"■^-ri'm.rn "^.™ 



± i li^rn rh 






+ ... 



+ .. 



, (3.80) 



where the operators M^ ,^ and E^^^ define the corresponding residues and ellipses denote the 
contribution of the k—th order poles {k < N — 2). 

The operators M^.Tn and E^,yj share all properties of the operator Q+(u,u). In particular, 
they commute with each other and are diagonalized by the eigenstate \E'p-^{q^g}(z). Replacing 
the operators Mm,m and E^,^ by their corresponding eigenvalues, Rm,rh{Q^Q) aiid E^jf^(q,q), 
respectively, one obtains from (|3.8CI| ) the pole expansion of the eigenvalues of the Q+— operator, 
Qq^q{u, u). The properties of the eigenvalues, Eq. ( |3.51| ), are translated into the following relations 



■± 



Rtrn{(l, Q) = (-1) e'^-' Rl,rn{-(1^ -^) ' ^m,,n(9, Q) = -^^,^(-g, -q) • 



(3.81) 



Similarly, substituting u = u^ + e and u = u^ + e into ( p.55| ) and matching Qq,q{u^ + e, u+ + e) 
into ( p.80| ), one arrives at 



r: 



El 






Ju;„.n r^N 



' C^^ [Q,Ms + m-l),t{s + m- l))X 



Nem,fh - i-T- In [Qq,q{i{^ + m - 1) + e, z(s + m - 1) + e)]" 



e=0 



(3.82) 



-Cs/[Ti'^B{m, 1 — 2s)B{m, 1 — 2s)] is expressed in terms of the Euler S— function, 

= iIj{1 — 2s + m) + iIj{1 — 2s + m) — iplm) — ip{m) — 



the constant C. is defined in (3.561) and e 



ijj{l - 2s) + i'{2s). We verify, using ( |3.82D and (|3.45|) , that -R^i(g, q) = e''""''' Cf in accordance 
with dpI). 

Remarkably enough, the spectrum of the model can be found from the pole expansion of the 
Q+— operator, Eq. ( |3.8CI| ). Indeed, according to ( p.81| ), the quasimomentum 9q^q can be expressed 
in terms of the residue functions Rfi{±q,±q). Substituting ( p.80|) into ( p.74|) , we obtain that 
the energy E]y{q, q) is related to the residue functions E^^{±q, ±q) 



EN{q, g) = - Re [E^^,{q, q) + E^^^{-q, -q)] = Re [E+i(g, q) + E+^^-q, -q 
Making use of ( p.82|) it is easy to see that this expression coincides with ( p.73|) . 



(3.83) 
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3.6. Asymptotic behavior 

As we have seen in the previous Section, the eigenvalues of the Q— operators, Q'q^q{u,u), are 
meromorphic functions of the spectral parameters, u = X—in/2 and u = A+m/2, with a (infinite) 
series of poles in A located parallel to the imaginary axis outside the strip | Im A| < 1/2. Let us 
find their asymptotic behavior at large A along the horizontal axis. Re A ^ oo for Im A = fixed. 
In this limit, we may neglect the imaginary part of A and choose the spectral parameter to be 
real, A = z/, with z/ -^ oo. 

We substitute u = —in/2 + u and u = in/2 + u into the relation ( |3.5|) and examine its both 
sides in the limit of large real u. Since the Q— operator commutes with the SL(2, C)— generators, 
Eq. ( |3.37| ), it can not depend on the momentum operator p and the same is true for its eigenvalues, 
Ql^q{u, u). Using this property we put p = in (|3.5| ). The reason for this particular choice of the 
momentum is that the state \E'p-=o,{g,g}(^) = ^{gg|(^) possesses additional symmetry properties - 
it diagonalizes the operators 5*0 and S^o, and at the same time it is annihilated by the operators 
p = iS- and p = iS- 

S. M/|°;-j(f) = , (So-h) v^;°^;_j(z) = O (3.84) 

and similar relations hold in the antiholomorphic sector. These properties imply that, firstly, 
\E'l -t(z) is translation invariant and, as a consequence, it depends on the differences of the 
coordinates and, secondly, it has a definite scaling dimension 

vl/(o)(^, + e, Zfc + e) = vl/(0)(z,, z^) , ^^'\\zk, Izk) = A'^-^^A^-^^"vI/(o)(^,, z,) . (3.85) 

Let us determine the asymptotic behavior of the operator Q+(m, u) on the space of functions 
^(^''(i') satisfying the conditions ( p.85| ). Using integral representation for the Q+— operator, 
Eq. (|3.6| ), with the kernel given by the second relation in ( |3.27D , one arrives at the following 
relation in the limit z/ ^ oo 



N 



[u,u 



|vl/W](i') ^=- u'^^'-^--^^ c{srs) U^z, - z.^i]'-'' (3.86) 

k=l 

N 



k=l 



with c(s, s) = [TTa{2 — 2s) (—1)"=]^. Examining the asymptotics of the integrand at large u, one 
finds that the dominant contribution comes from two different integration regions: 

(I): Wk = 0{u), (II): Wk-Wk+i = 0{l/u) (3.87) 

with A; = 1, .., A^, in which the //—dependence of the integrand cancels out in the product of the 
z— dependent factors. 

In the first region, at large Wk, one rescales the integration variables as Wk = Vk^-' and Wk = Vk^ 
and applies the identity 

[uyk-Zkriuyk-Zk+i]-''' ^ exp{-^{zk-Zk+l)/yk-^{zk-Zk+l)/yk) =e-^'''<^^''-''^^^^^^^ . (3.88) 

Taking into account the scaling properties of the functions '^^^\ Eq. ( p.85| ), we obtain from ( |3.86|) 

_(m,m)^(°)]j(z) -^°- j,h+-h-Nis+,) ^ cis,s)^'\z), (3.89) 
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where the subscript (I) refers to the integration region in (|3.871 ). Here, the notation was in- 



troduced for the function ^J depending on the differences of the coordinates Zk — Zk+i- It 
is obtained from the function '^^'^\z) through the transformation generated by the operator Si 
defined as 

¥^\z) = [Si¥^^]{z) (3.90) 

fc=i -^ "^ fc=i 

where Vy = (fyi...(PyMS{^-^^yi/N) is the integration measure on the space of translation 
invariant functions. According to (|3.90|) , the operator Si transforms \I/*^°^(z) into another trans- 
lation invariant function. Then, it follows from ( |3.89| ) that the contribution of the region (I) to 
the asymptotic behavior of the operator Q+(m,-u) on the space of the functions satisfying ( |3.85|) 
is given by 

[Q+(Z/ - m/2, Z/ + m/2)]i "=°° lyh+h-Nis+s) ^ ^^^^ -^ ^^ ^g^g^^ 

with the operator Si defined in ( |3.9(]| ). 



In the second region, Eq. ( p.87|) , all Wk approach the same point zq, but its position on the 



plane can be arbitrary. It is convenient to introduce the "center-of-mass" coordinate Zq and 
define the relative coordinates y^ 

Zq = jriwi + ■■■ + wn) , iJk = ^{wk - Zo) = 0{l''^) , (3.92) 



such that J2k Vk = 0- Changing the integration variables in (|3.86|) from {wi, ..., w^r) to {zq, yi, ..., y^) 
one applies the identity 

[zo + Vk/v - Zk]'"'[zo + yk-i/v - ZkY"" ~ exp(-i(|/fc - yk-i)/{zk - Zq)) (3.93) 

and uses the translation invariance of the function ^^^^ to find after some algebra 

[Q+{u,u)¥%^{z) ""^ z/i-'^+i-^-^(^+^")xc(s,5)5i?)(z). (3.94) 

Here, \l/jj (z) is translation invariant function related to '^^^^{z) through the transformation 
generated by the second operator Sn defined as 

5;?)(i') = [5„^(°)](f) (3.95) 

AT „ AT „ 

k=\ "^ fc=l "^ 

with the measure Vy defined in ( ^.9(J| ) . Then, the contribution of the region (II) to the asymptotic 
behavior of the Q+— operator can be expressed as 

[Q+(z/ - zn/2, V + zn/2)]jj "=°° ^^'h^-x~h-N{s+s) ^ ^^^^ -^ ^^^ _ ^3_gg) 



We notice that this relation looks similar to ( |3.91| ) and it can be obtained from the latter by 
replacing the conformal spin of the state, h ^ 1 — h and h ^ 1 — h, and substituting the operator 
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Si by its counterpart (|3.95| ). Moreover, comparing ( |3.90|) and ( p.95|) one finds that the kernels of 
the operators are related to each other as 



Si{z\y) 



N 

n 

fc=i 



N 



Zk " Zk+i 



l-2s 



Su{y\z) Yllvk 



Vk+i 



\2s-l 



(3.97) 



fc=i 



with the operator of mirror permutations M defined in ( |2.58| ). 

Finally, combining together two contributions, Eqs. (|3.91| ) and (|3.96|) , we obtain the following 
asymptotic behavior of the operator Q_(_ 



t+ 



{u — in/2, V + m/2) 



c(s, s) 



^/i+/i-Af(s+S) ^ _|_ ^l-/i+l-fe-Af(s+S) ^ 



(3.98) 



We would like to stress that this relation is valid on the space of functions satisfying the con- 
ditions ( |3.85|) .P| One finds the properties of the operators 5i/n by requiring that the asymptotic 
expressions for the Q— operators have to satisfy the relations established in Sect. 3.2. In partic- 
ular, it follows from ( 3.29 ) and ( |3.42 ) that these operators commute with each other and with 
the Q— operators 

[cSi, Q±(u, u)\ = [Su, Q±{u, u)] = [Si, Sii] = . (3.99) 

Let us compare (|3.98D with the general expression for the asymptotic behavior of the eigen- 



values of the operator Q+(z/ — m/2, v + in/2) that follows from (|3.55| ) 
g|+)_j(i/-m/2,i/ + m/2) "=' 



(+) I piWq^q/2 piN{ips-nns)/2 
{9.?}' 



Q\[i,\ Q-^l^li-Q-^i-" -■-'11 " y 



-2iNus 



(3.100) 



with s = {l + ns)/2 + iiys and the phase (ps defined in ( |3.55| ). Matching ( |3.98| ) into this expression 
we require that the z/— dependent part of the phase should not depend on the conformal spin of 
the state h = {1 + nh)/2 + iuh and h = 1 — h*. This condition is satisfied provided that the 
eigenvalues of the operators Syu have the same absolute value and differ only by a phase 



>I/II ^q,q 



\Zi,Z2,...,Zn 



T,.e'(^--±®--)v^S(^"'i,^-'2 



9:5^ 



zn 



(3.101) 



operator (|3.98|) have the following 



As a consequence, the corresponding eigenvalues of the Q 
asymptotic behavior 

Q^g^jiu - in/2, V + in/2) ''=°^ 2^^-^'^^+^'''^^ c{s, s) T^^^e^^^-' cos (8^,5 + 2uh\nu) (3.102) 



with Tqq, Qgg and flg^g being real. Its comparison with (|3.100|) yields (up to an additive 
(g, g)— independent correction to Qq^q coming from the normalization factors) 



Q 



q,q 



-Waa ■, 



Q 



(+) 

{g,g} 



~ v 



l-N 



COS (Og^q + 2z//, In z/) 



(3.103) 



We recall that Wq^q was defined in (|3.55|) as eigenvalue of the operator W. 

Since eigenvalues of the Q— operator do not depend on the momentum of the state p, the 
asymptotic behavior ( |3.102| ) holds for arbitrary eigenstate '^ j^^i^q^q^{z) . This is in distinction with 
( p.98|) that is valid only on the space of functions satisfying the conditions ( p.85| ). 



This implies, in particular, that the states belonging to this space have an infinite norm ( |2.9| ) and therefore 
one can not define the operators conjugated to Si/n with respect to the scalar product (2.9). 
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In this Section, we have constructed the Baxter Q— operator and have shown that the main 
physical observables of the system, hke Hamihonian, quasimomentum operator, transfer matrices 
etc, can be expressed in terms of a single operator Q_|_ («,«). Thus, the problem of finding the 
energy spectrum of the model is reduced to solving the second order finite-difference Baxter 
equations, Eqs. ( |3.3| ) and (p.4|), on the eigenvalues of this operator. The construction of the 
corresponding eigenstates \E'p>^{g^g}(ii, ..., ijv) will be the subject of the next Section. We would 
like to stress that solving the Baxter equations on the eigenvalues of the Q— operator, we have to 
impose the additional conditions on their solutions that follow from the analytical properties of 
the Q— operator, Eq. ( |3.8CI| ), and its asymptotic behavior at infinity, Eq. ( ^.1021 ). For a general 



solution to the Baxter equation, depending on the total set of the integrals of motion, {g, g}, 
these conditions can be satisfied provided that the values of {q, q} are quantized. The explicit 
form of the quantization conditions will be discussed in the forthcoming publication [|^ . 



4. Separation of Variables 

For the system of A^ = 2 particles, the eigenstates can be found exactly, Eq. ( |2.55|) , due to the 
SL{2, C) invariance of the Hamiltonian. For arbitrary A^, due to complete integrability of the 
model, the functions \E'p'^{g^g}(i') can be defined as simultaneous eigenstates of the total set of the 
integrals of motion, p and {qk,qk) with k = 2,...,N. The latter are given by the k—th order 
differential operators acting on holomorphic and antiholomorphic coordinates of the particles. 
Instead of going through the solution of the resulting system of the differential equations on 
'^p,{q,q}{z) we shall employ the method of the Separation of Variables (SoV) developed by Sklyanin 



11]. It allows to find the integral representation for the eigenstates of the model by going over 



to the representation of the separated coordinates x = (xi, ...,xn_i) 

'^f,{q,q}{z) = / dx fi{xi, ...,Xn~i)Uj^,xi,...,Xn_^{Zi, ..., Zn) {^{q^q}{Xl, . .. , XN-i))* , (4.1) 

where Uf^g{z) is the kernel of the unitary operator corresponding to this transformation 

Up,xi,...,XN_i{Zi, ..., Zjy) = {Zi, ...,Z]y\'p,Xi, ...,Xiy^i) , (4.2) 

and ^^qgy{xi, ...,XAr_i) is the eigenstate of the Hamiltonian in the SoV representation 



^{q,q}{Xl, -, XN-l)5ip - p') = (^p',{5,g}|p',fl,...,f7V-l) = / rf2;?7^',5-i,...,xjv-i(^) (%,{9,g-}(^)) • 

(4.3) 
Here, we introduced the standard notation for the bra and ket vectors on the quantum space of 
the system and defined the eigenstates in the different representations by projecting them out 
onto the corresponding basis of the states, (\E'p'_{g^q}|p',xi, ...,xn-i) and (\&p-^{g^g}|ii, ..., zn)- 



Remarkable feature of (^4.1| ) is that its substitution into ( p. 81) leads to the Schrodinger equation 
for ${gq}(x) that takes the form of multi-dimensional Baxter equation (see Eq. ( |4.19| ) below). Its 
solution is given by the product of the eigenvalues of the Baxter Q— operator with the spectral 
parameters equal to holomorphic and antiholomorphic components of the separated coordinates 

(<l>{g,g}(xi,...,X7v_i))* ~ Q{xi,xi) . . .Q{xn-i,xn-i) ■ (4.4) 



^°Here, it proves convenient to define the transformation $ ^ ^ to be anti-linear. 
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In this Section, we shall construct explicitly the SoV transformation, Eqs. ( [4.1|) -( ^^ , establish 
the quantization conditions on the possible values of the separated coordinates Xk and calculate 
the integration measure /i(xi, ...,xn-i) entering ([4.11). 



4.1. Representation of the Separated Variables 

Defining the representation of the Separated Variables for the noncompact SL(2, C) magnet, we 
follow closely the Sklyanin's approach [|lT|. This approach is based on the properties of the auxil- 



iary monodromy operator Tj\j{u) defined in ( p^.44| ) and, in particular, the off-diagonal component 
Bi\i{u) and its antiholomorphic counterpart. One finds from the Yang-Baxter equations on the 
monodromy operator that the operator Bm{u) satisfies the following relations [jl], ^ ^ |1T| 

[B^{u),Bn{v)] = [S-,Bm{u)] = 0, [So,Bjv{u)] = -Bn{u) (4.5) 

with Sa being the total S'L(2,C) holomorphic spin. 

We recall that, according to Eq. ( p. 45 ), Bn{u) is given by a polynomial in u of degree A^ — 1 



with the operator valued coefficients. Commutativity property ( [4.5|) implies that the operator 
coefficients commute with each other and, therefore, can be diagonalized simultaneously. The 
eigenvalues of the operator B]sr{u), being polynomials in u, can be parameterized by their zeros 
Xi, ...,xn-i- Denoting the corresponding eigenfunction as U^^xi^i, ..., z^) and taking into account 
( p.45|) , one gets 



Bn{u)U^^^{zi, ..., zn) =p{u- xi)...{u - xn^i)Uj^,x{^i, •••, zn) (4.6) 

and similarly in the antiholomorphic sector 

Bn{u)Uj;^s{zi, ..., Zn) = p{u - xi)...{u - xn-i)Uj},x{^i, •••, zn) ■ (4.7) 

Here, x = (xi, ..., xn-i), with Xk = {xk, Xk), denotes the zeros of the eigenvalues of the operators 
Bn{u) and Bn{u)- 

The common eigenfunctions of the operators B^iu) and B^iu) define the basis on the quan- 
tum space of the system, U^^^^Zi, ...,zn), which is parameterized by the momentum p and the 
set of zeros x . Using this basis, one may expand the eigenstates of the Hamiltonian, \l/p' {^ g} (i*) 
over the states U^^g{zi, ...,zn) to arrive at the expansion similar to ( [4.1| ) |^. As was shown by 



Sklyanin |TT| , the resulting expansion defines the representation of the Separated Variables. 

Since the eigenfunction in the separated coordinates, Eq. ( ^4.4| ), is symmetric under permuta- 
tion of any pair Xk and Xj, we impose the same condition on the solutions to (^^ ) and ([4.7|) 



'Jp^xx.--x,....Xj...Xf^^i\Z) — Up,xx...Xj...Xf,...xi^_i\Z) ■ 14. Oj 

The operators Bn{u) and Bn{u) are conjugated to each other with respect to the 5'L(2,C) 
scalar product, Eq. ( 2.46|) . Together with ( [4. 61) and ( |4.7|) this implies that the same relation holds 



between their eigenvalues leading to 

xl = Xk, k = l,2,...,N-l. (4.9) 



^^ Taking into account ( ^.46l ), one notices that this basis consists of the eigenstates of two mutually commuting 
hermitian operators, B{u) + B{u*) and i{B{u) — B{u*)), and therefore it is expected to be complete. 
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together with p* = p. This becomes the first constraint on the possible values of the separated 
coordinates. As a consequence of (|4.9| ), the solutions to ( ^^ and (|4.?1) corresponding to different 
Xk and p are orthogonal to each other with respect to the scalar product ( p.9| ) 

{x,p'\x,p) = j £zU^^^{Zi,...,ZN)\U^^^,{Zi,...,Z]^yj 

= i2nf6^'\p-p'){6{x-x') + ...}^f^, (4.10) 

where ellipses denote the sum of terms involving all permutations of the vectors inside the set 
X = (xi, ..., xn-i) that are needed to restore the symmetry property ( ^4.8| ). To give the meaning 
to the 6{x — x') one has to specify the possible values of the separated coordinates (see Eqs. (|4.35|) 
and (i;^) below). 

4.1.1. Sklyanin's operator zeros 



Following [11|, one can interpret parameters Xk in ([4.6|) as eigenvalues of certain operators Xk 



and represent the operator Bn{u) as 

Bn{u) = iS^{u — Xi) . . .{u — xn~i) ■ (4-11) 

One also defines similar operators Xk in the antiholomorphic sector, so that Bi^{xk) = Bj^ixk) = 
0. According to ( |4.5| ), the operator zeros x^, defined in this way, satisfy the commutation relations 



[xk.Xj] = [S^,Xk] = [So,Xk] = (4.12) 

and, therefore, they can be diagonalized simultaneously 

^k Ujj^xi^) = Xk U^^xi^) , iS^ Ujj^xi^) = p Ujj^xi^) . (4.13) 



To find the Schrodinger equation on the expansion coefficients ^^q^q}{xi, ...,xn-i) in ( [4.1|) we 
introduce the operators X^ 

X+ = ANixj ^ u) , Xr = DNixj ^ u) , (4.14) 

where Xj <-^ u stands for the substitution of the spectral parameter by the operator from the 
right [jll| (see Eq. (|4.20| ) below) .P^ The operators X- and their antiholomorphic counterparts 
X- , J = 1, ..., X — 1, have the following remarkable properties [|ll|] 

X^ Xj = {xj T iSjk) X^ , X^ Xj = {xj T i^jk) X^ (4.15) 

and 

^k U,^,xi^ = (xk ± ^s)''U.^^^^^^iz) , X± U.^iz) = (x, ± ^sfU^^^^^^g^{z) , (4.16) 



^^In [Q, similar operators, X^, were defined by substituting u = Xj in ( [4.14|) from the left. Both sets of the 
operators are related to each other as Xj^ = X Jliaifel^fc ~ % ^ i)/{ik — ij) and it is more convenient to our 
purposes to use the operators X^. 
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where the vectors e^ and e^ add a unit to the holomorphic and antiholomorphic components of 
the vector Xk = {xk,Xk), respectively, 

x±iek = {xi,...,Xk±iek,...,XN) , 4 = (1,0), (4.17) 

and similar for 5 ± i e^ with e^ = (0, 1). Then, the equation on ${g^g}(a;i, ...,xn-i) follows from 
the following operator identity 

{^p,{q,q}\tNiXk ^ U)\p',x) = tNiXk){^f;,{q,q}\p',x) 

= {xk + is)^ {'^p,{q,g}\p' , X + i ek) + (xk - is)^ {^jj,{g,q}\p',x - isk) , (4.18) 

where, in the first relation, we took into account that \E'p> {g^^} (i*) diagonalizes the auxiliary transfer 
matrix tj^{u) = A]\r{u) + Dj^{u) for arbitrary u and, in the second relation, we substituted 
tNixk •«—>«)= X^ + X^ and applied ( [4.161 ). 

Substituting ( ^l3|) into the last relation, we find that the expansion coefficients ${q,q}(5) have 
to satisfy the multi-dimensional Baxter equationj^ 

tAr(xfc)${q,q-}(5) = {xk + is)'^^{g^g}{x + iek) + {xk-is)'^^{q,q}{x-iek) (4.19) 

with A; = 1, 2, ..., A^ — 1 and t^lxk) being the eigenvalue of the auxiliary transfer matrix (|2.47|) , 
depending on the quantum number g^. Similar relation holds in the antiholomorphic sector. 
Thus, in the representation, defined by the functions U^xi^^ ^^^ Schrodinger equation on the 
states '^p,{q,g} (-?) becomes equivalent to the system of separated, one- dimensional finite difference 
equations on ${gq}(£). It remains unclear, however, what are the possible values of separated 
coordinates x and how does the integration measure fi{x) in ( [4.1| ) look like. Both questions will 
be answered in the next Section. 

One can find the recurrence relations on the integration measure yu(5) by taking into account 
the properties of the operator An{u) under conjugation, Eq. ( |2.46|) . Since An{u) is a polynomial 



in u of degree A^, it can be uniquely reconstructed from its values sd, u = Xk with k = 1, ..., A^ — 1, 
Eq. ( [4.14| ), and the asymptotic behavior at large u, Eq. ( p.45| ). 



N-l N-1 



AN{u) = {u + tSo + J2^k)lliu-x,)+X^J2Il^-IY- ^^-^^^ 

k j=l k=l j^k ^ ^ 

Writing similar expansion for An{u) and imposing the condition A]^{u) = [Aj^{u*)]\ one obtains 
the relation between the operators [X^Y and X^ 

iXfy = X^Yl \-lf\ (4.21) 

Here, we also included the relation between [Xj^y and Xj^ that comes from the analysis of the 
operator Dn{u). Combining together ( [4.21J ), ( [4.13| ) and ( |4.16| ), one obtains 



'±\trr /r?\ _ /^ _1_ „•„^Af TT ^'^ Xj±l 






Xk Xj 



^•^ Since x HSk does not satisfy ( [l.g| ), arri ving at this equation one assumes certain analyticity properties of 
the function <&{g,g} (x) outside the region ( |4.9| ) that wiU be discussed below. 
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Then, examining the identity 

Jd'z {x^U^^,,x'i^)yU,^,xi^ = j<fz [U.,^^,{z)y (X±) V^.,^(z) , (4.23) 



we calculate its both sides using ( [4.16|) , (|4.13|) and ( [4.21|) , and compare the coefficients in front of 



5{x — X ±iek) in the resulting expressions. In this way, we arrive at the following finite-difference 
equation on the integration measure 

^i{x ± iSk) ^ -r-r Xk -Xj±i 

This relation constraints the dependence of the measure on the holomorphic coordinates Xk- One 
gets similar relation in the antiholomorphic sector by replacing Xk ^ Xk- 

The solutions to the equations (|4.19| ) and (|4.24|) are defined up to multiplication by an ar- 



bitrary periodic function /(5), such that f{x ±iek) = f{x). The algebraic approach described 
in this Section, does not allow to fix this ambiguity and, therefore, the construction of the SoV 
transformation remains uncomplete. To avoid this obstacle, we shall construct in the next Section 
the exphcit expression for the basis functions U^xi^)- 

4.2. Construction of the SoV transformation 

The transition functions to the SoV representation, U^^{z), are simultaneous solutions to the 
system of equations ([4.6|)-( ^rT0D . Let us start with the conditions ( |4.6| ) and (|4.7|). Since Bn{u) 



is a polynomial of degree A^ — 1 in the spectral parameter u, its general expression at arbitrary 
u can be reconstructed from its special values at A^ — 1 distinct points that we choose as -u = x^. 
Then, Eqs. ( |4.6| ) and ( |4.7|) become equivalent the system of equations 

BN{xk) Ujjxizi, ..., zn) = BN{xk) Uj^xizi, ...,zn) =0, k = 1,2, ...,N - 1 . (4.25) 

One of the reasons for this particular choice of the spectral parameters is that similar equations 



have already appeared in our analysis of the Q— operator, Eqs. ( |3.19| ) and (|3.18|) . 
Indeed, using (|3.18|) and ( p.l9|) , we find that the function 

Uf,xi^) = / d^y ^i'i'Si(^i' •••' ^N I m, •••, yiv)^p,x2,...,xjv^i(y) 



\im[yi] ^' / rf\yj';^J(ii,... ,£7V 1^1,^2, •••,yW)^p,x2,...,x^_i(^) (4.26) 



satisfies the relations ( |4.25|) at A; = 1 for arbitrary weight function Z{y) and y = {y2, ...,yi^). It 



remains to show that the same relations hold for k > 2 and, in addition, to ensure the symmetry 
property ([4.8|). It is easy to see that the former requirement is equivalent to the latter property 



combined with Bn{xi) U^r(.{z) = Bn{xi) U^^{z) = 0. Thus, to satisfy the equations ([4.25|) for 
k > 2, it is enough to require that the integral in r.h.s. of ( ^.26| ) should be symmetric under 
permutations Xi ^^ Xk- This condition imposes constraints on the possible form of the weight 



function Z{y) in ( |4.2(j| ). Supposing that such weight function exists, one applies the last two 



relations in Eq. (|3.19|) to obtain 

AN{xk)U^^xi^) = {xk + is)^U^^£^i^^{z), 

DN{xk)U^^,xi^) = (xfc-^s)X,*-ie,(^), (4-27) 
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where the vector e^ was defined in ( [4.17|) . These relations coincide with the operator identities 

o. 

As a hint to finding the solution for ^p,x2,...,ijv-i5 ^^ notice a similarity between ( [4.26| ) and the 
definition of the kernel Xy^y.^u,u{S\ w), Eq. ( p.32| ). For v = Xi and v = Xi both expressions involve 
the same F— function but defined for different values of yi. Most importantly, the X— function 
is symmetric under permutation of the spectral parameters, Eq. ( p.35|) , and it is this property 
that we want to impose on ( [4.26| ). Matching the second relation in ( [4.26| ) into ( p.32| ) we choose 
the weight function as 



^p,X2,...,XM-l^ 



{a{s + ix2,s -1x2)) 



N-l 



dwA 



(l-s,l-s) 



{y2,---,yN\w3,...,WN)z, 



p,X3,...,XN-l 



.W . 



(4.28) 

Substituting this ansatz into ( |4.26|) one finds that the dependence of U^j^^z) on xi and X2 is 
factorized into the following integral 



{s,s) 



(l-s,l-s) 



^iV-l(^l)^7V-2' (^2) izi,...,ZN\w3,...,WN) 



(4.29) 



{s,s) 



(l-s,l-s). 



d yKuxA^li •••' ^N I Z/2, ..., l/Ar)A^2,S2 (1/2, •••, VN I W3, •••, Wat) , 

where subscript in the l.h.s. refers to the number of right arguments of the A— function. We would 
like to stress that the number of left arguments of the function Ivxl'^l^ is larger by 1 than the 
number of its right arguments. Therefore, in the convolution of two A— functions this difference 
increases to 2. Diagrammatical representation of ( |4.29| ) is shown in Fig. |1^. The symmetry of 
( [4.26|) under permutation of Xi and X2 comes from the following remarkable property of (|4.29| ) 



A(/4(^i)A!^:f" 



"■N- 



[X2] 



a{s + ixi, s — ixi^ 
a{s + ix2, s — 1x2] 



N-l 



A%'%{X2)A 



{l-s,l-s) 
Af-2 ' 



Xi, 



(4.30) 



which can be considered as a generalization of the similar property of the X— function, Eq. ( |3.35| ). 
Replacing (s, s) — *> (1 — s, 1 — s) in ( [4.30|) , one obtains another useful relation 



A 



{l-s,l--s) 
N-l ' 



(s,5) 



^l)A7V-2(^2 



a[s + ixi, s — ixi^ 
a(s + ix2, s — 1x2] 



l-N 



A 



(l-s,l-S)/-j X a(s,s) /-j X 
N-l [X2)^'^N-2\^^) 



(4.31) 



which is valid for arbitrary N, Xi and X2- 

The proof of ( ^.30| ) and ( ^.31| ) goes along the same lines as those for the X— function, 
Eq. ( p.35|) , and it is based on the permutation identities shown in Figs. || and |^. Namely, 
one inserts two auxiliary lines with the indices ±i{xi — X2) into the leftmost rhombus in Fig. |T^ 
and moves one of the lines to the right of the diagram, systematically applying the permutation 
identity. Fig. |^. In contrast with the previous case, the chain of rhombuses is not periodic and 
two auxiliary lines can not annihilate with each other and remain attached to the right- and 
leftmost rhombuses. Nevertheless, each of these lines disappears due to the identity shown in 



Fig. 13 



Applying the identities ( [4. 301) and (|4.31 



it becomes straightforward to write the expres- 
which is consistent with ( |4.26| ) and is symmetric under 
permutation of any pair of the separated coordinates, Eq. 



sion for the transition function U^^{z 



U, 



p,X\ 



Cn{X 



){pY-'^" fd'wNe'^^^-^- U^{z;wn) 



(4.32) 
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Figure 12: Diagrammatical representation of the convolution of two A— functions defined in 
Eqs. ( ^:29D . 



a(s+iv) a(s- iv 




1-s+iv 



a(s+iu) a(s-iu) 

1-s-iu 




Figure 13: Special case of the permutation identity shown in Fig. |^ and obtained by sending one 
of the external points to infinity. 



where 2(p ■ wn) = pwN + PWn, 



for even N, and 
U^{z;wn) 



(4.33) 






p,X\ 



in agreement with ( [4.13|) . 



l(fi)A!i:f-^"^(f2)AK(^3)...AS'~'''"'^(xiv~i)](^i,...,iW|^iv) (4.34) 

for odd N. Here, the convolution involves the product of {N — 1) functions Ki^^-^ixk) with 
alternating spins (s, s) and (1 — s, 1 — s). The following comments are in order. 

Since the A— functions are translation invariant, their convolution in ( [4.33|) and ([4.34|) depends 
on the differences zy, — wn and, as a consequence, (iS'_ —p)U^ 
The additional factor {p'^y^ 

Throughout this section we have tacitly assumed that A^'^^'|2j.(^i, •••, ^Af | ^2, •••, ^Af) is a well- 
defined function of the z— and y— vectors on the plane. As we have seen in Sect. 2.2.1, this 
condition constraints the possible values of the spectral parameters, Eq. ( p.21|) . Repeating the 
same analysis for the function A^^J^^.^, defined in (|3.18| ), we find that the separated coordinates 
have to satisfy the relation ( |2.21| ) with u and u replaced by Xk and Xk, respectively. Together 
with (|4.9|) this leads to the following quantization conditions on the separated coordinates 



in (|4.32| ) restores the scaling dimension of f/p-^-H 



Xk = J^k 



ink 



Xk 



ink 
^k + ^ 



(4.35) 



with Vk real and Uk integer. We would like to remind that it is for these values of the spectral 
parameters that the i?— matrix is a unitary operator, Eq. ( p.38|) and (|2.39| ). 



^"^ Since the separated variables are dimensionless (see Eq. ( 4.12| )), it follows from ( 4.10|) that the scaling dimen- 
sion of C/-.34 is equal to (A'' — 1). 
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Figure 14: Diagrammatical representation of the transition function, U^^{zi, ...,zn)- The mo- 
mentum p flows into the diagram through the top vertex wn. The indices ak and jSk parameterize 
the corresponding two-dimensional propagators (see Eq. ( |A.1| )) and are defined differently for even 
and odd k: a2fc = 1 - s - ix2k, /52fc = 1 - s + ix2k and a2k+i = s - ix2k+i, P2k+i = s + ix2k+i- 
Integration over the positions of (A^ + 1)(A^ — 2)/2 intermediate vertices is tacitly assumed. 



The normalization factor cn{x) in ([4.32|) ensures the symmetry property ([4.8|) of the expres- 
sion ( ^4.32| ). It compensates the additional factors in the r.h.s. of ( [4.3U| ) and ( 4.31| ) and is given 
for A^ > 3 by 



[(7V-l)/2] 



[iV/2-1] 



cn{x) 



Y[ [a{s + ix2k, s - ix2k)]^ '' JJ [a{s + ix2k+i, s - ix2k+i] 



(4.36) 



k=l 



k=l 



while C2{xi) = 1. Using the definition of the a— function, Eq. ( [A.3| ), together with the relations, 
s = 1 — s* and Xk = xl, one finds that the r.h.s. of ([4.36|) is a phase factor, |c7v(5)| = 1- Its 
explicit expression at A^ = 3, 4 looks as follows 



C3[Xl,X2) 
Ci{Xi,X2,X3) 



[a{s + ix2, s - ix2)] , 

[a{s + ix2, s - 1x2)] [a{s + ix^, s - ix^)] 



(4.37) 



The expressions (|4.33|) and (|4.34|) can be represented as the "Pyramide du Louvre" [^ 
diagram shown in Fig. 0. This diagram consists of {N — 1)— rows, in accordance with the total 
number of A— functions in Eqs. (|4.33| ) and (|4.34|) . The A;— th row of the pyramid is built from 
the lines carrying the indices ak and (3k that depend on the separated coordinate Xk and are 
defined in Fig. |14|. It is attached to the {k + 1)— th row through (A^ — A;)— vertices. At A^ = 2 the 
expression for the pyramid looks like 



(4.38) 



Up,a-M = \p\ I d'w2e''P-^' UsM^2) , 

U,,iz; W2) = [z, - «;2]^-"H^2 - «;2]^+*"^ 
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S-IX. / \ S+IX1 



^1 *v 




/ ^^ 


1-s+ix'^ \ 


W2 


/ 1-s-ix'^ 



Figure 15: Diagrammatical representation of the scalar product of two pyramids at A^ = 2. 

At iV = 3 the corresponding expression for [/^^^^^^{z; W3) is equal to the product of two (nonunique) 
star diagrams, which can be expressed in terms of 2-P'i— hypergeometric series. We do not present 
here its explicit form, since it is more convenient to our purposes to use the integral representation 



4.3. Integration measure 

The transition functions U-j^{z) defined in ( [4.32|) are the eigenstates of two mutually commuting 
hermitian operators, B{u)+B{u*) and i{B (u) — B {u*)) , Eqs. ( [4.6|) and (|4.7] ). As such, they should 
be orthogonal to each other with respect to the SL{2, C) scalar product (|2.9|) . The general form 
of their orthogonality condition, consistent with the symmetry properties ( ^.8| ), is given by ( |4.1CI| ) 



and it involves the integration measure fi{x). To find fj,{x) we substitute the transition function 
( [4.32| ) into the orthogonality condition (|4.10|) and match the result of integration into the r.h.s. 
of (gJ^). 



Let us evaluate the scalar product ( [4.1(]| ) using the diagrammatical representation of U^xi^ 
shown in Fig. 0. The conjugated function {U^j^{z))* is represented by the same pyramid dia- 
gram, in which the exponents a^ and j3k are replaced by the corresponding conjugated exponents 

ttfc -* ttfc = 1 - ttfc , (3k ^ (3k =1- (3k- (4.39) 

It becomes convenient to flip horizontally the conjugated pyramid diagram, so that the point W]^ 
will be located at the bottom of the diagram and the points Zk at the top. The scalar product 

(x.».|S'.»W = /<;=.[/^(f;».)(C/^.(£;.-W)* (4.40) 

is obtained by sewing the pyramid and its conjugated counterpart through the points Zi, ...,zn. 
The resulting Feynman diagram has the form of a big rhombus built out of [N — 1)^— elementary 
rhombuses (see Fig. |16] below) with the top and bottom vertices located at the points wjy and 
w'j^, respectively. Notice that the indices a^ and Pk in the upper and lower part of this rhombus 
depend on the separated coordinates x and x, respectively. 

Let us first calculate (|4.4CI|) at A^ = 2. Substituting ( [4.38|) into ( [4.40|) one obtains the rhombus 



diagram shown in Fig. [T^. Integration over zi and Z2 can be easily performed using the chain 
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relation (|A.6|) , leading to the expression independent on W2 and w'2, {xi,W2\x'i,w'2) ~ |w2 — 



W2|°. Its Fourier transformation with respect to W2 and ■W2 leads to a singular expression.P^ 
To identify the corresponding generalized function, one has to regularize the underlying two- 
dimensional Feynman integrals. The singularities appear due to the fact that the sum of four 
indices, corresponding to the sides of the rhombus, coincides with the space-dimension, dg = 2. 
This suggests to regularize the Feynman integrals by shifting the indices as 

(ofc, ak) -^ (a/c -e, ak-e) , {(3k, (3k) -^ {(3k -e, (3k- e) (4.41) 

with 6 real positive. The sum of four indices inside the regularized rhombus becomes 2 — Ae and 
the Feynman integrals are well-defined. Their calculation is based on the chain relation ( [A.6| ) 
and it leads to 

{xi, W2\x^,W2)£ = vr^[u72 — Wg]"^^ a{s — ixi — e,l — s + ix[ — e,l + i{xi — x[) + 2e) 

X a{s + ixi — e,l — s — ix[ — e,l — i{xi — x'l) + 26) , (4.42) 

where the subscript e in the l.h.s. indicates the regularization. Going over to the momentum 
representation and using ( [A.2| ), we calculate the Fourier transform of ( [4.42[ ) and obtain the 



expression involving the factor a{—Ae) = T{1 + Ae)/T{—Ae) that vanishes as e — > 0. To get a 
nonzero result, the smallness of a{—Ae) should be compensated by one of the factors in 



Examining (f4.42| ) one finds that this happens at xi = x[ since a(l + 2e) = r(— 2£:)/r(l + 2e) ~ 
1/e.Q Carefully examining the limit e —>■ one obtains 

lima(-4e, l + i{xi - x[) + 2e, l-i{xi -x[) + 26) = 27c6n.n'J{ui - u[) = 27t6{xi - x[) , (4.43) 



where xi = i^i — ini/2 and x[ = v[— in'i/2. Finally, combining together Eqs. (|4.42|) and ([4.43|) , 
we obtain the orthogonality condition at A^ = 2 



{x[J\xiJ) = \iuif / d2w2ci'w^e2™2P-='™2P'(f^,^2|x'i,w^)e = {2Tif5^^\p-p')5{xi 



TC 



4 



X, 



(4.44) 
Its comparison with ( [4.10| ) yields the expression for the measure at A^ = 2 

lj{xi) = 2/7r\ (4.45) 

The calculation of the scalar product at higher N goes along the same lines and is shown 



schematically in Fig. 0. Examining the rhombus diagram shown in Fig. p!6|a, one notices that 
all sewing vertices Z2, ...,zn-i, except the left- and rightmost ones, zi and 2; at, respectively, have 
four lines attached to them. At the same time, only two lines are joined at the points zi and 
zn- The integration over these two points can be performed using the chain relation, similar 
to the previous case, producing two vertical lines with the indices ±.i{xi — x'l) (see Fig. |l^b), 
which can be moved horizontally through the diagram in the direction of each other, thanks to 
the permutation identity shown in Fig. 0, until they meet and annihilate each other. In this way. 



one arrives at the diagram shown in Fig. 16c. It differs from Fig. 16a in that two sewing vertices 



^^This has to d o with the fact that ( 4.40 ) should be understood as a distribution, that is one has to integrate 
the r.h.s. of ( 4.4C ) with a smooth function of x and perform the integration over z afterwards. 

^^There exist other possibihties to get the factor ^ 1/e by putting, for instance, xi = i{l — s) for xi ^ —is, 



but the corresponding values of the separated coordinates do not satisfy the quantization condition (4.35) 
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(a) 



(b) 



(c) 



(d) 



(e) 



Figure 16: Diagrammatical calculation of the scalar product of two pyramids. The fat points 
indicate the vertices that can be integrated out using the chain relation shown in Fig. 



zi and zm disappeared and the spectral parameters xi and x[ were interchanged. We notice 
that in this diagram there are already four vertices which can be integrated out using the chain 



relation ( [A.6|) . Further simplification amounts to repetition of the steps that we just described. 
Each next step reduces the number of vertices in the diagram and interchanges the spectral 



parameters Xk 



x:- 



Continuing this procedure, one obtains the diagram shown in Fig. p!6|d, in 



which {N — 1) elementary rhombuses are aligned along the vertical axis and are attached to their 
neighbors through the common vertices. Each of these rhombuses corresponds to the Feynman 
diagram in Fig. |15| defining the scalar product at A^ = 2, Eq. (|4.42| ). The only difference with 
the previous case is that the separated variables parameterizing the A;— th rhombus from the top 
are given by XN-k and x^ instead of xi and x[. We expect that the chain of rhombuses should 
produce the contribution ~ Y[k=i^i^N-k — x').). To calculate it carefully, one has to regularize 



the Feynman integrals according to (|4.39|) . Applying ( [4.42 ), we replace each rhombus by a single 
line with the index —Ae (see Fig. ^Be) and integrate over their (A^ — 2)— common vertices using 
the chain relation ( |A.6|) . Going over to the limit e —>■ and taking into account the identity 
( [4. 431 ), we combine together different factors coming from ( [4.3(3| ) and ( [4. 421 ), perform the Fourier 



transformation ( [4.32| ) and finally obtain the following relation 



lim {x,p\x ,p') 



e^O 



(27r)^ 6^'\p - p') (5(2) (xi - x'^_i)5(2)(x2 - x'^_,) . . . 6^^\xn^, - x[) 



N~l 



X 



-vr^ Y\_ a{l + i{xk-Xj),l-i{xk- Xj)) . 



(4.46) 



j>k 



In arriving at this relation we integrated out, using the chain relation ( |A.6| ), all vertices of 
the pyramid except the top and bottom vertices. Each integrated vertex brought the factor 
~ 7ra(l + i{xk — x'jy-j)) and their product appeared in the second hne in ([4.46|) . 

Comparing ( [4.46| ) with the general expression for the scalar product, Eq. ( [4. 101 ), we notice 
that it is not symmetric under permutations of the separated variables x^ ^-> Xj. The reason 
for this is that performing the calculation of ( [4.46| ) we have tacitly assumed that xi is different 
from x[, X2,. 



x'j^_2, thus eliminating all terms in ([4.10|) except the one entering ( [4.46| ). Indeed 
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as we have seen in the case of the N = 2 scalar product, the hmit xi —>■ x'l requires the 
regularization of the exponents ( [4.41D . Although this regularization makes all Feynman integrals 
well defined, it transforms the unique Feynman diagrams (stars and triangles) into nonunique 
ones and, therefore, it does not allow to apply the permutation identity (see Fig. |). Since in 
our calculation we applied the permutation identity to all pairs of the separated coordinates 



{xi,x'j^_i^), k = 2, ...,N — 1, except {xi,x'j^_i) we have to assume that \xi 



X 



N-k 



7^0. 



Matching ( |4.46| ) into ( [4 .101 ) we obtain the following expression for the integration measure 

-7V2 N-l o^-N^ ^-1 



^{x) 



2ti- 



(AT-l)! 



J]^ \a{l + i{xk 






1-2 



2ti- 



j>k 



(N-l)\ 



n 

j>k 



\Xk 






(4.47) 



which is valid for A^ > 3. Here, \xk — x^p = {xk — Xj){xk — Xj) and the relations ( |A.3D and 



have been taking into account. Replacing Xk by their quantized values (|4.35|) one gets 



fi{x) 



2tt- 



-Af2 N-l 



{N-l)\ 



n 

j.k^i 

j>k 



[vk-i^jf + -{nk-njf 



(4.48) 



We verify that the obtained expression for the integration measure satisfies the relation ( [4.24| ). 

Having determined the integration measure in the SoV representation, we can use ( [4.10 ) to 
write the completeness condition 

Jd'pjdxf,{x) (U.^{z',,...,z'j,)y U.^^{zu...,z^) = {2nf6'^'\z-z'), (4.49) 

where the integration over quantized values of the separated variables, J dx, implies the summa- 
tion over integer n^ and integration over continuous Uk defined in ( |4.35| ) 



N-l 



dx 



n E 



k=l 



\nk-- 



dvk 



(4.50) 



Since the eigenstates \E'p'^{q^g}(i') are orthogonal to each other with respect to the S'L(2,C) 
scalar product ( p.9|) , the same should be true for the functions <I>{g^g}(xi, ...,X7v-i) in the SoV 
representation. To obtain the corresponding orthogonality condition one substitutes ( [4.1| ) into 
]9D and uses Eq. ( [4.10|) . In this way one arrives at the Plancherel formula 



\^P:{'?.<?}I^P'.{9'.<?'}/ 



{2Tiy5{p-p'){^ 



{q,q}\'^{q',q'}l 



) = {27rr6{p-p')6{q-q'), 



where the notation was introduced for the norm of states in the SoV representation 

i^{q,q}\^{q',q'}) = / dx fx{x) (<l>{g,g}(f i, ..., f ^_i)) *${g,,g-/} (f i, ..., f 7V-1) = 5{q - q') . 



(4.51) 



(4.52) 



Here, ${g,g}(xi, ..., xn-i) is a completely symmetric function of the separated variables Xk, satisfy- 
ing the mult i- dimensional Baxter equations, Eq. ( |4.19|) , in the holomorphic and antiholomorphic 
sectors. As we will show in Sect. 4.4, the solution of these equations is factorized into the product 
of the eigenvalues of the Q— operators that have been studied in Sect. 3. Therefore, substitution 
of ( [4.4|) into ( [4.52 ) yields the orthogonality condition for the solutions of the Baxter equation. 
Its explicit form depends on the normalization factor present in the r.h.s. of ( [4.4[ ) and it will be 
established in the next Section. 
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4.4. Relation to the Baxter Q— operator 

The construction of the SoV representation, performed in Sect. 4.2, is similar in many respects 
to that of the Baxter Q— operators in Sect. 3. In addition, the eigenfunctions in the SoV rep- 
resentation, ${gg}(a;i, ...,a;Ar_i), and the eigenvalues of the Baxter Q— operators satisfy similar 
finite-difference Baxter equations. This suggests p^ that there should exist the relation between 



two different objects - the transition function to the SoV representation, U^r[.{z), defined in 
( [4.32|) , and the Baxter operators Q±{u,u) given by Eqs. ( p.27|) and (|3.28| ). In this Section, we 



shall establish such relation (see Eq. ( 4.54| ) below) and use it to prove 



To start with, let us consider the product of A^ — 1 Baxter operators Q_ (m, u) 

= j Sy... d'y'Qi-X{z\ y')... Qi~l,,,^_,{y\w) , (4.53) 

where y = {yi,...,yN) and the kernel Qlj^^l{zi, ...,Z]s[\wi, ...,wn) was defined in ( p.28| ). Using 
the diagrammatical representation of the Q^xi,xi{z\y'), the left diagram in Fig. ^ we associate 
( [4.53| ) with the Feynman diagram, in which (A^ — 1)— periodic chains of the unique triangles (one 



for each Q_) are glued together through their common vertices as shown in the left diagram in 
Fig. |T^. Comparing this diagram with the one shown in Fig. |T^, we notice their similarity - 



the lowest row defining the i*— dependence look alike, as well as their difference - the number 
of vertices is higher in the former diagram; the latter diagram does not have periodic boundary 
conditions along each row. In addition, we notice that Q^{z\ wi, ...,wn_i,wn) depends on the 
{N — 1)— additional arbitrary vectors Wi, ...,wj\f^i as compared to U^{z; wjy), defined in (|4.33|) 



and ( [4.34| ). As we will see in a moment, in order to obtain the transition function Ug.{z] wn) from 



( [4.53|) , these vectors have to be chosen in a particular way - one has to put wi = W2 = ■■■ = wn- 



and take the limit wi — >■ oo afterwards 

Up,xi^) = Cj; Jim [w;i]2(^-i)(i-^) / (fwN e'^^^+'P^^ Q^{z\ w,, ..., ^i, ^;v) ■ (4.54) 



uii— >oo 



Here, Cf is the normalization factor depending only on the momentum p. 

The proof of ( [4.54| ) is based on the reduction formulae obtained in Sect. 3.3.4 and it goes as 
follows. The dependence of Qgi,...,£j^_j^{z\ w) on wi, ...,W]y_i is carried by the rightmost Q— kernel 
in ( [4. 531) , or equivalently by the A^ unique triangles in the upper row in Fig. |1^. At Wi = W2 = 



... = Wiv_i, (A^— 2) of these triangles shrink into 5— functions (see Eq. ( |A.9| )), leading to significant 
simplification of the corresponding kernel. Indeed, systematically applying the reduction formula 
( |3.58| ) one finds that Ql~^_^,SM-iiy\^) ~ "^(^i " yi)H^i - y2)--S{wi - yN-2)- Substituting 
this expression into ([4.53|) , one finds that the adjacent kernel Q^XM-2,XM-2^y" I V) ^^ evaluated at 
yi = y2 = ■■■ = yN-2 and the same reduction formulae can be applied again reducing the number 
of unique triangles in the row next to the upper one by A^ — 3. Repeating this procedure one 
finds that the substitution wi = W2 = .■■ = wn-i in the upper, (A^ — 1)— th row of the diagram, 
creates an avalanche of simplifications throughout the whole diagram, in which [k — 1) unique 
triangles disappear in the A;— th row {k = 2,3,...,N — 1) as shown in Fig. |l^. The resulting 
diagram still depends on the arbitrary vector wi, which defines the position of the vertex of the 
right- and leftmost unique triangles in each row. In the limit Wi —>■ oo, these triangles scale as 
a power of [wi], so that each row of the diagram contributes the factor [wi]~^*^^~'^^ Combining 



together the factors coming from (A^ — 1) rows of the diagram, one finds that the r.h.s. of (|4.54|) 
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Figure 17: Reduction of the product of {N — 1) Q— operators, Eq. (|4.53|) , into the pyramid 
diagram. The A;— th row of triangles in the left diagram is given by Fig. |^ for u = x^- The 
triangles shown by the dash-dotted lines collapse into the 5— functions at Wi = W2 = ••• = wn-i 
due to the reduction formulae. The triangles shown by the dashed lines are replaced by a power 

of [wi] as ifJi ^ oo. 



has a finite limit as wi ^ oo. The corresponding Feynman diagram has the form shown in the 
r.h.s. of Fig. 0. Still, it differs from the pyramid diagram. Fig. |14|, by additional horizontal lines 
with the indices (2s — 1). 



To demonstrate the equivalence between two different diagrams. Figs. [T^ and |14|, we apply 
the same trick as was used in Sect. 3.2 in establishing the commutativity condition for the 
Baxter Q— operators (see Fig. H). Namely, we start with the diagram in Fig. |l^ and insert two 
horizontal lines with the opposite indices, (2s — 1) and (1 — 2s), into all rhombuses with the 
indices {ai, Pi,a2, P2), in the lower part of the diagram. Applying the permutation identity. 
Fig. 1^, we move {N — 2) auxiliary lines with the indices (2s — 1) vertically to the upper part of 
the diagram until they either reach the boundary of the pyramid and disappear in virtue of the 
permutation identity. Fig. ^, or end up in the rhombus with the indices {aN-2, I3n~2, C(n-i, Pn-i) 
just below the top vertex wn- In the latter case, one gets rid of the auxiliary line by applying 
the star-triangle relation to the unique triangle with the indices {a^-i, Pn-i, 2s — 1). Up to the 
additional line with the index 2(1 — s) attached to the top of the pyramid, the resulting diagram 
looks like the original pyramid with one row less. The latter has been replaced by the lowest 
row of the diagram in Fig. |1^. Applying the same transformations to the reduced pyramid with 
(A^ — 2)— rows, we reconstruct one more layer of the diagram Fig. 17. We continue this procedure 
iteratively A^ — 3 times until we arrive at the reduced pyramid containing a single rhombus. 
Inserting, as before, two lines with the indices (1 — 2s) and (2s — 1) into the rhombus, we get 
rid of the latter line by applying the uniqueness relation to the triangle containing the top of 
the pyramid. The resulting diagram coincides with the right diagram in Fig. |1^, in which the 
chain of the additional lines with the index 2(1 — s) is attached to its top vertex wn- The total 
number of lines in this chain, [(A^ — l)/2], is equal to the number of rhombuses in the original 
pyramid, in which two opposite vertices are located on the vertical axis going through the top 
vertex wn. Going over to the momentum representation, the whole chain is reduced to a power of 



the momenta ((— 1)""']9^ "^^p^ 2s^[(Af i)/2] ^^^^ ^^^ recover the relation ([4.54|) . Carefully collecting 
all factors and comparing the expression for the r.h.s. of ( [4. 541 ) with Eq. ( f4.32| ), we calculate the 
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normalization constant in ( [4 .541) as 



C.= ((_l)n.pl-2.-l-2s-y(A^-l)/2] [^^(2(1 - ,))]-(A^-l)(A^-2)/2 _ (4,55) 



It comes about from applying the reduction formulae (|3.58| ) to (A^ — 1)(A^ — 2)/2 collapsing unique 
triangles in Fig. ^ and compensates the difference of the scaling dimensions in the both sides of 
the relation ( [4 .541 ). 



The relation (|4.54|) between U^^{z) and the Baxter Q— operator can be rewritten in the 
operator form by introducing the special state \uj^^^^), which depends on two arbitrary vectors, 
Zq and p, and is defined in the coordinate representation as follows [^ 



N-l 



{z\u,,^,^) = c^ [zo]'(^-')('-^) ^iP^N+^v-ZM -Q 5(4 _ ^^) . (4.56) 



fc=i 



Using this state, one finds from ( [4.54| ) 



Up,x{^) = Ji"^ {z\Q^{xi,xi) . . .Q-{xn-i,xn-i)\uJzo,p) ■ (4-57) 



^Q— >00 



Many properties of the SoV representation follow from this relation. In particular, the sym- 
metry property ( [4.8|) comes from the commutativity of the Q— operators. Substituting (|4.57|) 
into ( ^4.3|) , we take into account that the eigenstates '^^^{q^q]{z) diagonalize the operator 



Ui,M 



^{q^q}{x)5{p-p') = (%{g,g}|Q_(Xi,2;i) . . . Q^{x N-l, Xn-i) \^zo=oo,p') 



N~l 



c^S^^\p-p')l[Q'q;^{xk,Xk), (4.58) 



fc=i 



where the normalization factor is defined as c$ = J cPp (\E'p,{g,g}| 1^^0=00,^' ) and Q^q^^{u,u) stands 
for the eigenvalue of the operator Q_(m, u). This relation replaces the ansatz ([4.4|) and it provides 
the expression for ${gg}(a;) that automatically satisfies the multi-dimensional Baxter equation 
( [4.19| ). Applying the conjugation relations ( p.40|) , one expresses the wave function in the sep- 



arated coordinates, ${g^g}(x) in terms of the eigenvalues Qg^q = Ql^q{u,u) (up to irrelevant 
normalization factor c<j,[a(l — 2s)/7r]^*^^~^)) as 



N-l 



${g,g}(:r) = Yl Hs + ixk, s - ixk)]'^ {Qq,q{xk, Xk))* ■ (4.59) 



fc=i 



We recall that the possible values of the separated coordinates, Xk = Vk— ink/ 2 and Xk = xl, are 
parameterized by integer Uk and real z/^, Eq. ( 4.35|) . The explicit form of ${g^g}(:r) at A^ = 2 can 



be found from ( |3.63|) . 

The properties of ${q^g}(:r) can be established using the results for the functions Qq^q{u,u) 
obtained in Sect. 3. The analytical properties of ^^q^qj{x) on the complex z/jt— plane follow 
from Eq. ( p.80| ) and Fig. |l^. In particular, ^^q^qj{x) does have any singularities on the real 



^^Notice that this state is not normaUzable with respect to the SL{2,C) scalar product (2^) and, therefore, 
it does not befong to the quantum space of the system. Keeping this in mind, we shah use the same scalar 
product notation for the convolution of the wave function ujgg^p{z) with an arbitrary test function, {^/\ujgg^p) — 
J (Pz{"i'{z))*ujSo,piz), assuming that the integral is convergent. 
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z/A;— axis and, applying ( |3.102| ), one finds its the asymptotic behavior at large z/^ and fixed n^ as 
${g^g}(5) = (9(z/^~^*^'*'^"'°'), leading to /i(a;)|${g^^}(5)p ~ ^/^l- These properties are in accord 
with the fact that ${qg}(:r) should be normalizable with respect to the scalar product ( [4.52| ). 
Substituting ( [4.5S ) into (|4.1| ), we obtain the integral representation of the eigenstates of the 



system, \E'p'|qq}(i'). As shown in the Appendix C, this representation can be rewritten, thanks 
to remarkable properties of the transition function U^ri,{z), in a concise form [0 



N-l 



^P,{9,g}(^) = n Qg,g(^k,Xk) ^p{z) , (4.60) 

where 



k=l 



N 



Qpiz) = / d'zo e^'P-'o -Q (^^ _ zor^^Zk - zo)"'' (4.61) 



fc=i 



is the so-called pseudovacuum state and Qq,q{xk,Xk) stands for the eigenvalue of the Baxter 
operator Q+ (u, u) with the spectral parameters substituted by the operator zeros. The repre- 
sentation ( [4. 601) is well-known in the theory of the 5'L(2,R) spin magnets. There, the relevant 
solutions to the Baxter equation are given by polynomials in the spectral parameters and ( 4.60|) 



becomes equivalent to the highest- weight representation for the eigenstates in the ABA approach. 
Eq. ([4.60|) provides a generalization of the above representation to the noncompact SL{2, C) mag- 
nets, for which the ABA is not applicable and the Baxter equations have nonpolynomial solutions 

iOT Qg^g{U,U). 

As a byproduct of (^39|), we substitute the relation ([4.59|) into the scalar product ( [4.52|) and 
obtain the orthogonality condition on the space of solutions to the Baxter equation 

dxn{x) YlQg',g'{xk,Xk) {Qg,g{xk,Xk)Y = 5{q - q') , (4.62) 

■^ k=i 



where the integration measure j dx ^{x) was defined in ([4.47|) and ( [4.50|) . 



5. Conclusions 

In this paper, we have studied completely integrable quantum mechanical model of N interacting 
spinning particles in two-dimensional space that has previously appeared in high-energy QCD 
as describing the properties of the A^— gluon compound states in the multi-color limit. Applying 
the quantum inverse scattering method, we identified this model as a generalization of one- 
dimensional homogenous XXX Heisenberg spin magnet to infinite-dimensional principal series 
representation of the 5'L(2,C) group. 

The model has a number of properties in common with two-dimensional conformal field 



theories |^. Introducing holomorphic and antiholomorphic coordinates on a two-dimensional 
plane one finds that the the kernels of the /?— matrix and the Baxter Q— operator are factorized 
into the product of two functions depending separately on the holomorphic and antiholomorphic 
coordinates, while the Hamiltonian is given by the sum of two commuting operators acting in the 
different sectors. In spite of this, the Hamiltonian dynamics of the model can not be reduced to 
independent dynamics in the holomorphic and antiholomorphic sectors. The interaction between 
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two sectors is ensured by the condition for the kernels of the operators and their eigenstates 
to be well-defined functions on a two-dimensional plane. This condition leads to the additional 
constraints on the properties of the model. In particular, it leads to the quantization of the 
integrals of motion of the system, the separated coordinates as well as the spectral parameters of 
the i?— matrix and the Baxter Q— operator. Using the global SL{2, C) invariance of the model, 
one can represent the eigenstates of the A^— particle Hamiltonian, "ifq^q^zi —zq, Z2 — Z0, ..., z^ — zo), 
as correlators of A^ + 1 quasi-primary fields in a two-dimensional conformal field theory. The 
latter are given by the sum of products of the holomorphic and antiholomorphic conformal blocks. 
This structure is general enough and, as was shown in Sect. 3.3.5, it also holds for the eigenvalues 
of the Baxter Q— operator. This, in turn, suggests that the Baxter Q— operators constructed in 
this paper should be related to similar operators defined in the conformal field theory P^j . 

We have argued in this paper that all integrability properties of the model are encoded in the 
properties of the Q— operator, satisfying the Baxter relations. Namely, the Hamiltonian as well as 
transfer matrices are expressed in terms of a single operator Q{u, u) and, therefore, their spectrum 
is determined by the eigenvalues of the Q— operator. In addition, going over through Sklyanin's 
construction of the Separated Variables we have demonstrated that the unitary transformation to 
these variables is also determined by the same operator. This leads to the integral representation 
for the eigenstates of the model with the wave function in the separated coordinates given by the 
product of the eigenvalues of the Q— operator. 

Thus, the spectral problem for the noncompact SL(2, C) spin magnet turns out to be equiva- 
lent to the problem of finding the eigenvalues of the Baxter Q— operator. The latter are defined as 
solutions to the Baxter equations, (|3.3| ) and (|3^ ), supplemented by additional conditions on their 
analytical properties and asymptotic behavior at infinity that have been established in Sects. 3.5 
and 3.6, respectively. The solution satisfying these conditions was found in the simplest case of 
the system of A^ = 2 particles. Constructing solutions to the Baxter equation for A^ > 3, we 
find that in order to satisfy the above conditions the values of the integrals of motion have to be 
constrained. The explicit form of the corresponding quantization conditions is known only for 



the polynomial solutions of the Baxter equation pT|, 30, 22|. Their generalization to arbitrary, 



nonpolynomial solutions will be presented elsewhere p4 |. 
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Note added 

When this paper was written, we learned about the recent work by H. J. de Vega and L.N. Lipatov, 
in which similar results were obtained - the unitary transformation to the separated variables 
was constructed, using different approach, for the SL{2, C) magnet of spin (s = 1, s = 0) with 
N = 2 and A^ = 3 sites; the explicit solution to the Baxter equation was found at A^ = 2 and its 
generalization to higher A^ was discussed. 
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A Appendix: Uniqueness relations 

Analyzing the noncompact spin magnet we encounter various two-dimensional integrals over the 
positions of particles which we represent in terms of the Feynman diagrams. Such interpretation 
turns out to be very useful as it allows to apply powerful methods of calculation the Feynman 
integrals well-known in perturbation QCD. More precisely, our calculations are heavily based on 
the "method of uniqueness" |l^, |16|. This method has been originally developed for the calcula- 
tion of d— dimensional Feynman integral and we shall use its simplified version at rf = 2. In this 
Appendix we collect all necessary formulae - the chain and the star-triangle relations, which form 
the basis of the method of uniqueness. Remarkable feature of the method is that subsequently 
applying these two relations one is able to obtain analytical expressions for complicated Feynman 
integrals without doing any integration. 

Throughout the paper we use the following compact notation for the "propagator" 

1 1 (^-w)"-" (;z-w)"-" (-1)°-" ,, , 



[z-w]'^ (z - w)°'{z - w)^ |2;-w|2° Iz-wp" [w - z]"" 

with a — a = n, and represent it graphically by the arrow directed from w to z and with the 
index a attached to it as shown in Fig. |1|. Going over to the Fourier transformation we define 
the propagator in the momentum representation 

Here, the notation was introduced for the function 

, , r(l — a) , , r(l — a) / n \ / \ / n\ / A ON 

a[a) = a{a) = , a{a, (J, ...) = a{a)a{P)... . (A.3) 

1 [a) 1 [a) 

It has the following properties 

a(l -\- a) 1 - f—i^a-a 

a(a)a(l-a) = l, A^ = __ , aia) = {-ir-^a) = \/ . , (A.4) 

aia) aa all — a) 



where in the last relation a — a has to be integer. Taking the limit of (|A.2| ) for a = a = ie as 



£ — i> we obtain the following convenient representation for the delta function 

5(2)(u;) = lim^— ^. (A.5) 

The method of uniqueness is based on the following two relations shown diagrammatically in 



Fig. ig and p 



Chain relation: 



where 7 = 2 — a — /? and 7 = 2 — 0; — /?. Ata + /5 = Q; + /? = 2 this relation can be further 



simplified using (|A.5|) as 



d^w- — = 7r^a(a,2-a)(5(2)(x-z). (A.7) 
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Star-triangle relation: 



7ra(a,/?,7) 



z - wY [x - wY [y - wX< [x - zY'^ [z - yy-P [y - x]i-° 



(A. 



with a + /3 + 7 = a + /3 + 7 = 2. 



a 



Z =71 a(a,p,2-a-P) t 



Figure 18: Tlie chain relation. 



a ,, 



= 7ia(a,|3,y) 1-P 




1-a 



a+|3+Y=a+|3+Y=2 



Figure 19: The uniqueness star-triangle relation. 



The l.h.s. of ( |A.8| ) is described by the "star" diagram in which three propagators are glued 
together at the same point w. The r.h.s. of ( [A.8| ) correspond to the "triangle" diagram. We 
shall call the star and triangle diagrams unique if, firstly, the sums of the holomorphic indices 
attached to three propagators in these diagrams coincides with the sum of the antiholomorphic 
indices and, secondly, it is equal to 2 and 1, respectively. Then, the star-triangle relation allows 
to replace the unique star diagram by unique triangle and vice versa as shown in Fig. ^. 

In the limit y ^ x, when the end-points of the unique star diagram approach each other, or 
equivalently the line in the unique triangle diagram with the index 1 — a shrinks into a point, 
one finds from (|A.7|) 



(fw- 



z — w]°' [x — w]i^ [y — wp 



'^='^7r^aia,2-a)5^^\z-x] 



(A.9) 



To prove ( |A.6D it is sufficient to take the Fourier transform of the both sides of the relation 
w.r.t. X and apply ( A. 2 ). The relation (|A.8| ) follows from ( |A.6| ). To see this one replaces the 
coordinates in (|A.6| ) as w -^ l/{w — y), x — > l/{x — y), z — > l/{z — y) and similarly for the 



antiholomorphic coordinates. After a simple calculation one reproduces Eq. ([A.8|) . 
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B Appendix: Relation between Q— operator and transfer 
matrices 



It is well-known in the theory of compact spin magnets that the transfer matrices of high spin 
are related to those of the lowest spin through the fusion relations []I8|, ^. Moreover, solving the 
hierarchy of the fusion relations one can express the transfer matrices of arbitrary spin in terms 
of the Baxter Q— operator ^^ |52|, ^. In this Appendix, we shall generalize these relations to 
the noncompact SL{2, C) spin chains. 

Our starting point is the relation between the transfer matrices and the product of Q— operator, 
Eq. ( |3.66| ). Using this relation we can establish some general properties of the transfer matrices. 

Since the SL{2, C) representations of the spins (s, s) and (1 — s, 1 — s) are unitary equivalent, 
we expect that the corresponding transfer matrices should be related to each other. Indeed, one 
finds using ( |3.66|) together with ( |3.67|) that the transfer matrices T^^°''^°\u, u) and T^^~*°'^~*o)(m, u) 
differ by a c-valued coefficient function 



T(so,so)(y^^) 



^{1-50,1-80)1 



M, U] 






a{s — 1 + So + iu, s — So — iu) 
a{s — So + iu, s — 1 + So — iu) 



N 



(B.i: 



Here, we should not worry about ordering of the transfer matrices due their commutativity. 

Making use of ( p.4(]| ), we obtain the expression for the transfer matrix ( 3.66| ) entirely in terms 
of the Q+— operator in two different forms 



¥("»'""») (u,u) 










= pJ;»''"«)(m,m) X [Q+{u* ~ tso,u* - iso)]^ Q+{u + tso,u + tso) 


(B.2 


= ^^°'''\u, u) X [Q+(n* - z(l - So), u* - t{l - sq))]^ Q+{u + z(l - sq), u + z(1 - 


- ^-0)) 


Here, the normalization factors are given by 




p^^''''\u,u) = 


\ -^^2i(u-u) "-s + 4:^s a{s-l + So + iu, s - So - iu)' 
TT^ a{s — So + iu,s + So — iu) 


N 


(B.3 


p<;-^"«)(«,«) = 


\i l)^^("-^)'''+/'''a(so s + tu,-s -So m)l 


N 







It is straightforward to verify that (p.2| ) satisfies the intertwining relation (p.lj ) and the normal- 
ization factors p^^°'^°\ p{so,so) a^]^^ p{so,so) g^^e related to each other as 



p^:"''\u,u) 



p(l-so,l-s-o)(^^ ^) M)^^^ M)/p(!-^°'l-^"°)(M, U) 



(B.4) 



Using (p.2| ) we can find the relation between the transfer matrix and its conjugate counterpart. 
Conjugating the both sides of the first relation in (|B.2|) we can either match it into itself using 



the relations (|3.49|) , or into the second relation in ( p.2|) . In this way, we obtain two equivalent 
representations 



[T^'°'''\u,u)] 



o-A]t 






J«o,so), 



—U*, — U" 



T^'°'"'\u* + i,u* + i) X 



(p(^°'^"°)Kw))^ 



Q 



p^Q°''°\u*+l,U*+l) 



(B.5) 
(B.6) 
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Finally, let us show that the transfer matrices of different spin satisfy nonlinear fusion rela- 
tions. To this end, one considers the product of two transfer matrices of different spins as well 
as different values of the spectral parameters and expresses it in terms of the Q— operators using 
( P.2|) . The resulting expressions involve four Q— operators including two conjugated operators. 
Using their commutativity, we may regroup them into two pairs of Q"''Q— operators and apply 



the relations ( [B.2|) to convert each pair into the transfer matrix. In this way, we arrive at the 



following fusion relation 

= f (^o-i-,so-ia)^^ ^^^^ + g) f('"+i"'^o+i^){u -g,u-g), (B.7) 

where T^*")('u) = T^'^°''^°\u,u)/p^^°''^°\u,u). Another fusion relation follows from the Baxter 
equation (|3.3| ). Multiplying its both sides by conjugated Q— operator and applying ( [B.2| ) we find 

tN{u) f("0'""o)(M,M) = (M + is)^f("°+^'^"«)L-^,MJ + {u - tsf T^''-'^'''^ (u + ^,UJ (B.8) 

and similar relation in the antiholomorphic sector 

iN{u) f^''''''\u,u) = (u + is)^f(^°'^"o+5)U,u--J + (u - is)^ f^^"'"""-^) U, u + - j . (B.9) 

There is the following important difference between these relations and ([B.7| ). According to 
the definition of the transfer matrix T^^°'^°\u,u), Eq. (|2.40|) , the spins (sq, sq) parameterize the 



auxiliary space V^^°'^°^ and for the SL{2, C) representation of the principal series they satisfy 
the conditions ( |2.12|) . We notice that for the transfer martices in the r.h.s. of (|B.8|) and ( p.9|) 
these conditions are not satisfied and, therefore, they do not belong to the same family of the 
transfer matrices of the SL{2, C) principal series. At the same time, appropriately choosing 
the parameters {g, g) and (cr, a), one finds from (|B.7| ) the fusion relations between the transfer 
matrices within the same family. 

C Appendix: Bethe Ansatz representation of the eigen- 
states 

In this appendix we prove the relation ( [4.60| ) which establishes the correspondence between 
( [4.1|) and the representation of the eigenstates in the Algebraic Bethe Anstatz |jl|, |^, |^. The 
ABA representation is based on the existence of the special pseudovacuum state fi(i'), which is 
annihilated by the spin operators of all particles (|2.2| ) 

N 

fi(^1 = n '>^''''k'' , 4'^ Q{z) = Si'^ Q{z) = . (C.l) 

fc=i 



The remarkable feature of this state is |jl], 0, ^ that it brings the Lax operators ( |2.15| ) to a 



upper-triangle form and, as a consequence, diagonalizes the auxiliary transfer matrices ( p.47| ) 

tN{u) n{z - zo) = [{u + isf + {u- isf] n{z - zo) (C.2) 
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with zq being arbitrary vector reflecting the 5*^(2, C) invariance of tN{u). The operator tAr(-u) 
satisfies similar relation. Then, the complete integrability of the model implies that Q{z — zq) 
diagonalizes the Baxter Q— operators and, as a consequence, the Hamiltonian of the model. This 
does not mean, however, that the pseudovacuum state can be identified as one of the eigenstates of 
the model. Examining its 5'L(2,C) transformation properties, Eqs. ( |2.1CI| ) and ( p.4| ), one verifies 
that Q{z — Zq) does not belong to the quantum space of the system (see footnote to Eq. (|4.56| )). 
Nevertheless, one can use its properties to establish different relations including Eq. ( |4.6CI| ). 
Going over to the momentum representation, we define the state 

npiz,,...,zr,)=c^^ Jdho e^'P^-'^^niz-zo) (C.3) 

with the normalization factor cO = (— l)"''p^*~^p^**~^/(7r'^a(2(l — s))) chosen for the later conve- 
nience. Let us examine its convolution with the transition function defined in (f4.57|) 



/ 



d^z{nj;{z)yU^,^{z) = hm {nj;\Q_{xi,xi) . . .Q^{xn-i,xn-i) \^zo,p') (C.4) 



zo=oo 



with the state co'?o=oo,p'(^) defined in ( ^.56| ). Then, substituting ( |C.3| ) into ( p.4| ) and calculating 



the convolution J (Pz{fl{z — zq))* Q^~^(i'|w), we represent the kernel of the Q— operator by the 
right diagram in Fig. |^. The integral over the intermediate points Zk takes the form ( [A.?] ) and 
it gives rise to the 5— functions which put the centers of the star diagrams at the same point Zq. 
The resulting diagram is given by the product of propagators connecting zq with the points z^ 
and it coincides, up to prefactor, with Q^. In this way, one gets 

{Qj^\Q^{u,u) = [iTa{2 — 2s,s + iu,s — iu)] {Q,^\ = A{u,u){Q^\ . (C.5) 

As it was expected, the pseudovacuum state diagonalizes the Baxter operator. One can verify 
that its eigenvalue satisfies the Baxter equation (|3.3|) with t^iu) replaced by (|C.2|) . 

Substituting ( |C5| ) into ( p.4| ) one replaces the Q— operators by their corresponding eigenvalues 
and reduces the r.h.s. of (|C.4|) to the convolution of the pseudovacuum state with the function 



LU 



zo,P' 



defined in ( ^.56| ). A simple calculation leads to 



Jim (fip-lt^^o^p-/) = Cjj6'-^\p - p) . (C.6) 



zo=oo 



Here, the normalization factor Cp- is given by ( [4.55 ). Finally, one calculates ( C4 ) as 



„ N-l 

/ d'z{npiz}rU.,,^{z)=c,^6^^\p-p')l[A{xk,Xk). (C.7) 

"^ k=i 

One can arrive at the same relation using the diagrammatical representation of the l.h.s. as the 
pyramid diagram. Fig. |14|, with the additional lines carrying the indices 2(1 — s) attached to its 
end-points Zk and joined at the same point ho- 
using the completeness condition, Eq. ( [4.49| ), we find from ( |C.7| ) 



s- 



/N-l 
dxix{x)U^^^^{z)ll{A{xk,Xk)T = {2n)-'' n^iz) . (C.8) 
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Here, the normalization factor in the l.h.s. compensates the difference in the scahng dimensions 
of f/p.^ and fip-. Let us consider an arbitrary operator f{xk,Xk), depending on the Sklyanin's 
operator zeros defined in Sect. 4.1.1, and apply it to the both sides of the last relation. Since 
f/p.^ diagonalizes the operators zeros, Eq. (|4.13| ), this operator can be replaced in the l.h.s. by its 



eigenvalue f{xk^Xk). Then, choosing f{xk^Xk) as the wave function in the separated coordinates 
${g,g}(£), Eq. ( [4.59|) , we obtain from ( [4.1|) (up to an overall normalization factor depending on 
the momentum p) 

/N-l N-1 

dXn{x)Uj..^{z)Y[ {A{Xk,Xk)yQq,q{Xk,Xk) = JJ Qq^q{Xk,Xk)flf{z) . 
fc=l fc=l 

(C.9) 

Here, Qq,q{xk,Xk) stands for the eigenvalue of the Baxter operator with the spectral parameters 
substituted by the operator zeros. 

The relation (|C9| ) generalizes the well-known highest- weight representation of the eigenstates 
of compact spin chain magnets in the Algebraic Bethe Ansatz 

|vI/fB^) = 5(Ai)...i?(A0|fip) (CIO) 



with B{u) being the off-diagonal element of the monodromy matrix, Eq. ( |2.44]) , and Ai,...,A, 



satisfying the Bethe equations. Using (|4.11|) , one can show that two representations, Eqs. ( P.9|) 



and ( ClOp , are equivalent if the eigenvalues of the Baxter Q— operators are given by polynomials 



of degree / in the spectral parameter u. Going beyond the class of polynomial solutions to the 
Baxter equation, one finds 0, ^ that ( |C.1CI|) is not applicable while the representation ( p.9|) 
remains valid. 
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